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ABSTRACT 


This report is a preliminary analysis of a ee 60ط‎ ۴67172066176167 
-. being developed at the Lawrence Radiation Laboratory, Berkeley. The 
spectrometer will be capable of presenting neutron spectra to an upper 
Et of 14.0 MeV. The spectrometer employs a EE gas mixture in 
& proportional counter, and the reactions 3He(n, p)T and EN d)D are 
utilized to obtain the incident neutron energies. Recoils from He, کا‎ 
and the stabilizing gas (CO2) are electronically removed from the spectra 
by taking advantage of the shorter rise times Of አርሽ ይጋ 6 ٠٣ 
analytic analysis of the rise time discrimination pulse rejection, the wall 
effects, and the spectrometer efficiency is presented. 

The computer program for the removal of the above corrections 18 
lengthy, and although an outline of the methods of machine computation 
15 presented, this ات‎ is not carried to fruition. Instead a computer 
program and an "inversion matrix" are developed (based on the analytic 
analysis M rented) by means of a "random probability analysis." The 
"inversion matrix" should permit the correction of a raw neutron spectrum 


output from a multichannel analyzer to a final neutron spectrum by 


-vil- 


time discrimination,‏ ا )5( ۲٦۹۶۰13 single operation for:‏ ۶< تحت 
(b) wall effects, and (c) spectrometer efficiency.‏ 

Appendices include a table of comparative physical and operational 
characteristics of all He proportional counters and ionization chambers. 
constructed, and reported in the technical literature to date (Fall 1965). 
The complete electronic circuitry necessary to construct the rise time | 
discriminator units, and the computer program ''CATMAN" developed 
for the random probability analysis are also presented as separate 


appendices. i 


eles 
Lo INIRODUCTION 
A. Background and Purpose 

This report concerning am iat oved 3He neutron SDeECirOmeuer eis 
a preliminary report concerning the theory, design, and construction 
of a ?He neutron spectrometer by the Health Physics Department of the 
Lawrence Radiation Laboratory. The actual operation and performance 
of the spectrometer will be the subject of a later report by the Health 
Physics Department and will be authored by Mr. Wai-Kit Quon. 

The detection of neutrons and especially the measurement of their 
energy have from the first attempts presented a more difficult problem 
Phan the other common nuclear particles. This difficulty is because the 
mea! methods depend on effects resulting from charges on the particles, 
. Most methods of neutron detection depend upon the scattering of particles 
after à collision with the neutron, and the measurement of the energy 
deposited by the scattered charged particles. The drawback to any 
method that is 09 to measure neutron spectra which is based on 
scattering lies in the conversion of the definite energy of the incident 
memenon to a continuum of energies of the recoiling particle. To avoid 
this continuum of energies, a neutron induced nuclear transformation 
was sought in which the energy release would correspond unambiguously 
to the neutron energy. Investigations of several potential reactions, 
and the elimination of those that had sharp resonances, low lying energy 


states in the residual nucleus, and large positive--or any magnitude 


negative O values, led to the selection of the following reaction: 


1 Í 


NECS —?H, E 


0 2 Hy + 0. 764 Me V. 


In addition, the 3He(n, p)T cross section is ey large and varies 
smoothly with energy, and the ability to discriminate against y rays 
enhances this selection for the task of neutron spectroscopy. Now, if 
the charged particles produced by the neutron interaction were to be 
completely stopped in the active volume of a proportional counter, the 
ionization caused by the loss of particle energy would be equal to the 
energies of the charged particles. This energy would then appear as 
free electrons and positive ions within the active volume. If an elec- 
trical field gradient is present, the electrons would migrate to the 
anode, and the charge deposited by the electrons would be proportional 
to the energy, E, of the incoming neutron. 

Several experimenters have examined the neutron interactions with 
the 3 He nucleus, and have constructed neutron detectors, and toga 
limited extent, neutron ዥች ም ማጮ a ae utilizing proportional counters 
and ionization chambers. The work of the majority of the previous 
experimenters has, in general, concerned itself with the neutron energy 
spectra of only afew MeV, although Sayres reports success with neutron 
energies up to 8.1 Mev!0 and the measurement of a cross section at 
17, 5 MeV vus using monoenergetic neutrons. 


In the majority of the previous work accomplished, two problems 


seem to inhibit the usefulness of a ?He neutron spectrometer. 


First, an ambiguity arises from the fact that ` He has a relatively 
large cross section for elastic scattering. The elastic 31e (n, . 
cross section is equal to the nonelastic cross section at . 072 Me v, but 
increases rapidly with increasing energy. At 2.2 MeV the els s R ros 
section is 4.6 times as large as the inelastic cross section. The ratio 
decreases slightly with further increasing energy, but the elastic cross 
Section is at least double the inelastic cross section from 0.25 MeV to 
14 MeV. The probability of an undesired recoil e is thus at least 
twice the probability of a desired inelastic collision. Unfortunately, the 
۶66011 events, in addition to being more probable, also produce a 
continuum of energy pulses from 3/4 E, to zero depending upon the angle 
Boat the سک‎ neutrons strike a ?He nucleus. This prevents utilizing 
any simple subtraction process in the determination of the initial neutron 
Spectra. 

second, the possibility of a nuclear event occurring in an orienta- 
tion such that the charged particles produced by the event intersect the 
walls of the containing vessel before expending their full energy was 
recognized, This phenomena gives rise to a count that is of less energy 
than a count that goes full track in the gas medium. ?He nonelastic 
collisions with high energy neutrons produce tracks that are of consider- 
able length in relation to the dimensions of the proportional counter, 7 
the wall correction effects increase significantly with increasing neutron 


energies. 


Solutions to the above two problems were attempted in various 
os 
manners. Brown constructed a ''wall-less'' counter composed 
of a peripheral ring of 16 small counters surrounding a central count- 
ing wire and an anticoincidence circuit to remove the wall effects. He 
then developed a computer program to unfold the He recoil events. 
lO 
Brown's work was extended by Wang who developed a computer programy 
to unfold the "wall effects" in addition to the "recoil events." Sayres 
10 ማቸ ae M 9۳ 

and Coppola developed a ''risetime discriminator" circuitry toe 
.removes the recoil events (and a fraction of the desired events) by taking 
advantage of the shorter track length of a 2e recoil in ء٦:‎ m 
the larger track lengths of the inelastic particles produced by neutron 
3 3 3,4 | 
interaction with à "He nucleus. Batchelor uses an approximate method 
of computation to remove the wall effects which appears to be useful at 
low energies, and short range of track in relationship to the counter 
dimensions. 

The purpose of this paper is to attempt to combine the removal of 


10 with the unfolding 


Ue ?He(n, n)?He recoil ceiicetdevcloped by Say ses 
oí the wall effect developed by ው and Wang. He It is further the 
purpose of this paper to report attempts to extend the energy range of 
the ?He proportional counter to 14. 0 MeV (which appears to be about a 
practical upper limit for these devices) and to develop a means whereby 


the spectrometer can be utilized with some degree of 4 x efficiency. ` The 


spectrometer will be used by the Health Physics Department, Lawrence 


Radiation Laboratory, Berkeley for radiation surveys--hence the 
requirement for 4x sensitivity. 
B. Description of the SECT Gele 

The experimental apparatus for this ?He 5066 0 ٦ ا‎ ete 
o a proportional counter, electronic components, and a computer 
Pogram. 
GS Proportional counter 

The proportional counter is a cylindrical tube with an active length 
of 15 inches, a 4-inch outside diameter, and a wall thickness of 0.035 | 
inch. The diameter of the anode wire is 0.003 inch. The proportional 
counter is designed with independently adjustable field tubes. Gas 
filling is 10 atmospheres of krypton, 2 atmospheres of high purity 3 He 
(less than 10-100 +۰٣٦ 17 and 1/2% CO» as a stabilizing gas. The 
tube was specially constructed by the Texas Nuclear Corporation for 


this spectrometer. A line diagram showing the tube and pertinent 


dimensions are shown in Fig. l. 
oS in, stainless steel tunin 
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A discussion of the design parameters and the design.variables 
J included in Appendix A. 
2. Electronic components 

The electronic components are designed to be able to distingui | 
between a pulse generated from a ?He recoil" and a He — — event 
To understand the operation of the electronics (a detailed description 
oí which appears in Appendix A), it is first necessary to under stand 
what the electronics are to accomplish. Basically, the jene of tram 
of a He recoil (or any heavy particle recoil), is dependent upon the 
mass of the incident particle, the mass and charge of the recoil par ٠ 
the energy of the incident particle, the scattering angle, and the composi- 
tion of the absorbing medium. (An analysis of the collision dynamics 
involved are included in Appendix B.) In general it can be stated a 
for equal recoil energies, a heavier particle, or a particle of the same 
mass but a higher atomic number, will traverse a shorter distance sS 
the same stopping medium than a particle of less mass. In a propor- 
tional counter, a track's radial component from the center wire determine 
the length 2 time that it takes to collect the charge, i.e. the pulse rige 
time. Since the track length of the recoils is always less than the track 
lengths of the disintegration particles (which are the desired eve ፲ቲ 


is possible to select a rise time that will effectively exclude all eec 


“The arrangement also serves to eliminate recoils from the krypton gas 
and the Cep 


events, yet retain a certain portion ot the 46 ٠ 7 ٤ +١٤١٠٣ 
have a longer rise time because oí their longer ር ንር. 
Unfortunately, the rise time criicrion also eliminates many ''good 
pulses'' whose tracks are orientcdu with respect to the center. wire in 
such a way that their rise times are within tme discrimination enitceuce 
mae electronic circuitry developed in Appendix Ais designed to 
accomplish this pulse rise time discrimination. 
Eu Computer program 

The computer program is developed to accomplish the "unfolding" 
of the wall effects and to restore those desired pulses e were | 
electronically eliminated by rise time discrimination. The wall effects 
computer program follows the method outlined by Wang, Ho but not 
ESenucal. Differences are primarily in the computation of the probability 
function P(E— Ei) which is the probability per unit energy that the 
reaction products of energy E will deposit energy E. inside the sensitive 
volume due to the wall effect. In addition, a new probability function is 
analytically developed to replace mathematically the valid pulses which 
were ቿን ው: due to rise time discrimination. An inversion matrix 
is proposed which will combine the two probability functions, and enable 
a complete reconstruction of the neutron spectra, above a minimum 
ያ that is determined by our rise time discrimator setting. ine 
neutron spectrum below this minimum energy cutoff level is unfortunately 


lost in the process of recoil elimination. 


یا ۔ 


| Unfortunately, the computer program analytically developed in 
(sss Cr euros so much computer time that it was not consider 
warranted to carry it to completion. Instead, a random probability 
analysis was carried out, and the inversion matrix was developed on 
the basis of these random probability calculations. An outline of the 


random probability method is included in Section II- E and Appendix E. 


E 


I. DEVELOPMENT OF THE °He SPECTROMETER 
A. Helium-3 and Neutron Interactions 
A neutron can interact with a ?He nucleus in five possible ways. 


These reactions which are energy dependent are as follows: 


3 


| He(n, p)T for neutrons of all energies 


a Eie for neutrons of all energies 
3. E d)D for neutrons of energies > 4.36 MeV 
4. Eo pn)D for neutrons of energies= 7.32 MeV 
5. Es 2 NE) EST 62067612 ጊስ ራሪሰሪ ንን ሺ seam 1. 

The cross sections for these reactions are reasonably well known 
in the energy range 0.001 eV to 14.0 MINES with the exception of the 
?He(n, pn)D for which no data are available. Itis assumed that tu 
exceptional case has à small cross section, even though it is energeti- 
cally possible. A comparative chart of the other reaction cross sections 
is included as Fig. 2. In addition to the em pn)D reaction, the 
` He (n, pan) ^H . reaction will not be considered in the determination of 
either wall effects or rise-time discrimination. The reasons for this 
Omission are: 
a. At the maximum energy expected to be encountered in this proportional 
counter (i.e. 14. 0 Me V), T (n, 2n) cross section is only 15% of the total 
non-elastic cross siete ፡ 
b. The neutron induced disintegration of the M nucleus involves four 


separate product particles--only two of which are charged. Four particle 


E 


dynamics are complicated, and the energy imparted to the various 
particles may span a wide range of variables. Note, however, that the 
threshold for this reaction is 10. 5 MeV and the kinetic energy of the 
charged particles produced by even a 14. 0 MeV incoming neutron cannot 
exceed 3.7 MeV. Since some of the energy of the incoming neutron will 
undoubtedly be carried away as kinetic energy of the two product neutrons, 
itis highly probable that the sum of the charged-particle kinetic energies 
will be less than 3. 7 MeV. Pulses from a 14.0 MeV neutron (n, 2n) 
reaction would thus appear on a multichannel analyzer at less than 3.7 
ቴ ስ. ፒ ሮሮ to Fig, 2, it is to be noted that the cross اہ‎ of the 
> He(n, p)T reaction is fairly constant between about 0.4 and 2.5 MeV, and 
the cross Section at these energies is at least 14 times as great as the 
(n, 2n) reaction at 14.0 MeV. Distortions of the 7 spectra would 
thus appear above 10.3 MeV at essentially 0% at threshold, up toa 
maximum of 15% at 14.0 MeV. This distortion would also manifest 
itself as a lower energy pulse, (<3. 7 MeV) but the allocation of reaction 
particle kinetic conse would result in an energy distribution that would 
not be unique,. but rather span from 0 to 3. 7 MeV. Fora conti 
neutron spectra this "lower energy distortion'' is not consider ed Ses 
due to the more probable (n, p) reaction at these lower energies which 
would essentially masx these degenerate pulses from higher energies. 

If we were to consider only ው interactions within the active 


volume of our container (i.e. no multiple.scattering events), and to 
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Energy Dependent ?He Neutron Cross Sections 


ZE 


likewise consider no neutron disintegration interactions with the materiali 
in the container walls, 2 would find em (after recoils are eliminated 

by rise-time discrimination and wall effects mathematically resolved) 
three reactions which would still tend to disturb our spectrum (irrespectiv 
of all corrections). These reactions therefore limit the maximum 
efficiency of our spectrometer. The three interactions are: 

| 1. The disturbance of the energy spectrum above 10.3 MeV and below 

የ ል. ስን ን ፒፓ ትምም ቢክ 

2. The disturbances of the energy spectrum above 7. 32 Me V anqa 
6.68 MeV from the ር ቢከ بر‎ reaction (if. this reaction does in fact occu 
3. The disturbance of the u spectrum at all energies due to the 
neutron induced krypton disintegrations. £ 

The cumulative error introduced by neglecting these reactions is 
unknown. 

The remainder of this paper is to be divided into four parts. The 
plan of "unfolding" the raw spectra to develop a complete neutron spectrun 
is as follows: 

Development of a probability function to replace the pulses eliminated 
by 'rise-time discrimination. ' (II-B) 


Development of a probability function to remove the wall and end 


This problem is discussed in detail in Appendix A, although no attempt 
1s made to correct the spectrum. 


شی 


effects. (II-C) 

Combination of the two correction factors to reconstruct the neutron 
spectra. (II-D) 

Reasons for and development of a random probability inversion 
matrix in lieu of the above II-D.  (II- E) 

11 IL should be mentioned here (ረ. CRE oI I dE کت‎ 
mathematical expressions for wall and end effects is primarily a restate- 
Mgeat Of the work of Dr. W. R. ን and ۷٦ Wang?!®, although 
some modifications have been made to the analysis of the later author. 

A number of their drawings, graphs, and formulas are repeated in this 
text for continuity and to facilitate A ` They will be designated 
Pena ‘dagger’ (+) if. attributed to Dr. Brown, and by a "double dagger’ 


(r1)if developed by Mr. Wang. 


E Development of a Probability Function to Replace Pulses 
Milmuinated by “Rise=time Diserimumarion| 


‘Asa result of the arbitrary orientations of the disintegration particles 
of the us. p)T and the ee d)D reactions, some tracks will be posi- . 
tioned such that the time difference between the arrival of the first 
electron and the last electron to the central counting wire will be quite 
short. This situation will occur when the radial component of the track 
lengths is small, or when the event takes place near the central portion 
of the counter where the electrical field gradient is high. In both of these 


orientations the count produced by a ''valid event" will be rejected by the 


' £74 


“As an approximation, it can be shown that by relating a epecifie particle 


Edu 


eleectromic circuitry as though it were a recoil pulse. It was consideraa 
that if an analytical expression could be developed for determining the 
time difference between the arrival oí the first and the last electron to 
the central counting wire, the ''rise time" of the pulses could be deter- 
mined mathematically. A probability function could then be developed 
based on the registered counts that would be representative of those good 
I و‎ ን 0خ 81۲۰۲667۴ کت‎ To develop this probability 
function, it is necessary however, to make four assumptions: 

(1.) The time required for the charged particles to deposit their 
energy in the gas medium is instantaneous. ٦ 

(2.) A nuclear event can occur with equal probability within any unit 
volume in the sensitive region of the proportional counter. 

(3.) The mobility of the electrons in the gas mixture is directly propor: 
tional to the field gradient over the limited range of the electrical field 
gradients in which we are interested. 


(4.) The electron mobility in the gas mixture (2 atmospheres Me 


l0 atmospheres Kr, and 0.5% CO5) can be computed by assuming that the 


energy E; = 3 mv;" ; solving for v,; and then TOURS an average velocity. 
+` ዖ > . ኣሪ 1 ] 
to track length (X.(E)) and solving for time t; = (WO) | I 


that the time required for a particle to deposit its a is about 3 nano- 
seconds for a proton (longest track length, highest velocity) of 10.0 MeV 
energy, and 0.6 nanoseconds for a 10.0 MeV 2He recoil particle (sho 
track length, smallest velocity). Rise times of interest to our counter | 
| 


Gircultry are in excess of 3 microseconds. 


NS 


electrons behave as ions, and that the mixture mobility can bevcomputed 
as Can the "mixture mobilities Of ONS. 

Assumptions (3) and (4) will be di^ issed in greater detail in the 
following c 
Eu iectron migration time from any point 7 ۷۶) ۶۶۷٦ 
metry. In a cylindrical proportional counter, tne electrical field 
gradient can be determined at any point r from the center of the counter 


by the relationship: ዣ 
ET) eee 


= 


H 


where E(r) is the field gradient at any radial distance r from the 


anode. 
7 = applied voltage 
r. = radius of the cathode 
nor Se radius oi the anode 


(This formula is developed in Appendix D.) 

Hor Our E. dimensions we can plot the field gradient at any, 
radius from the center wire by using the appropriate applied voltage and 
counter dimensions. Instead of directly plotting the field gradient as 
ET as of radius however, let us divide both sides of Eq. 1 by the 
total pressure of the counter filling gas in cm of mercury and plot this 
value instead. The reasons for this will become evident shortly. This 
plot can be seen in Fig. 3 for representative counter voltages. 


An electron produced by ionization at a radial position r will drift 


slowly towards the anode under the influence of the electrical field gradient. 


(volts /cin* cm-! Hg) 


fun 
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Pioo: 


Eis 


—ra =3.81x 1073 em re = 4.99 om 
=4. 
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Variation of Field Gradient/Filling Pressure (E/P) 
for Representative Applied Voltage 





ን 


As can be seen in Fig. ጋ, the electrical field grad CTE ን ስ ንም ٠ 
Beyond about 1 cm from the cathode varies very slow) see c R 
electron passes this point in its inward radial drift, the field gradient 
sharply increases. If we assume a direct proportionality between the 
drift speed and the -E/P ratio, we can obtain an expression in cylindrical 


E metry for the time it requires for an electron to reach the cathode, i.e: 


L 
E = eee (2) 
m(E/P) 

a 

TE n 
where m, the electron mobility is the slope of the speed vs E/P curve. 

PO 50 3 ۱ کو‎ 

Measurements of drift speed for ?He and Kr with 0. 595 CO, indicate 


that within the limited range of E/P which is the most important in deter- 
mining the drift time of electrons for cylindrical proportional counters, 
that this assumption of linearity does not appear to be unreasonable. 
Substituting the expression for field strength in the case of cylindrical 
coaxial MEE. we obtain the electron drift time from any radial 


position r. 


3 m 
tz E c . Ln 2c rods 
Vo II Vo Ta 
E ۲ a 
پت‎ Eel 
ra 
ete 
1 2 م۸‎ Te zl 
t= = kn — |> (3) 
pay ra |1ሪ | 4 
o e 
Gs) Vo La 


Z 
Emce -፦' is very small, the time.necessary to travel from the surface of 


mae anode to the geometrical center is negligible. 


-18- š 


The development of Eq. 2 and 3 follows the logic developed by 
Harling. 19 

We still must determine however, some means of combining ٣ 
mobilities of the electrons in a ''gas mixture" since data is available 
only for mixtures of individual noble gases and moderating gases. A 
theoretical development of this "gas mixture" mobility proble roki 
presented by Staub, *! however the necessary gas constants are not 
available to compute the final result. We must therefore make use of 
assumption (a) and assume that as a first approximation, the electron 
mobilities benave as do ionic mobilities. 

While it is realized that no simple theory can account for the electron 
diffusion, attachment, recombination, and agitation velocity, as the 
electron moves through a changing field gradient, it is considered that 
this approximation will cause at most an error of a constant factor. The 
electron mobilities € thus assumed to follow the same physical laws 
as do ions--although it is realized that the situation is considerably more 
complex for the fast moving electrons. 

Staub’s method of combining ion mobilitie 2 is as follows: 


pnmo EM x ° EN 
ሥ12 Po A P ጆዶ2 


where 41192 the ion mobility in the gas mixture 


the ion mobility in the component gases‏ ڈو رر 
partial pressures of the component gases‏ =2 را 
P = Pi t P5 - total gas pressure.‏ 


sus 


Substituting 1n (the eléctron mobility) tex p we have: 


v 


1 Pug ee 1 


= ——  . — +—— + — 4 
Inmix Ptot mie Fot MKr . (4) 





Using the sources explained in the footnote, " we find that for our tube 


0.3 cm/p-sec: per r 


ርን ema Hg) 
p 
i ያ پر رص‎ ርር መ : Hg) 
ےت‎ = 2.01 atmosphere 
Piz} = 10.05 atmospheres 
Substitution of the above values into Eg. € elds ٦ - 


TOND 


em/u-sec: (em = ena E Utilizing this value inm Com TE the appropriate 
counter dimensions we can bns migration time for an electron roma 
radial point in the counter. This has me accomplished in Fig. 4.for 
Ae sentative counter operating voltages. 

2, Rise time. Rise time of the pulse has been defined as the time 
Eimeerence between the arrival of the first electron at the central counting 


wire and the arrival of the last electron. By use of assumptions (c) and (d) 


m 


20 سس سس  "‏ چپ 
“Tne mobility factor for Kr was taken from the data of English and Hanna.‏ 


ES mobility factor for the helium is taken from the data of Friedes and 
Chrien, ll which is plotted in reference 19 (Harling). The reason for this 
use of Friedes and Chrien's data rather than Harling's is that Harling, 

(1) made his measurements using tHe with only small amounts oí ?He, and 
(2) he used substantial percentages of CO5. Although drift velocity of 
electrons in 2He appears to be rather insensitive to moderating gases, the 
Eoo content of Friedes and Chrien was considerably less than that utilized 
by Harling. 


Electron migration time (psec) 
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Fig. 4. Variation of Electron Migration Time from Various 
Radial Positions for Representative Voltages 
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we have been able to determine migration time, and have formulated 
an expression for the time for an electron to transit any radial distance. 


Thus the rise time (t...) can be expressed as 


|! ፲ጋ ደሬ 2 2 
t stj-t s =-ጨ ያደ ٦ - 5 
r 2 1 2 Ms D 5 3 (2) 


Mote that this is identical to the time that it takes an electron at 
position r. (furthest distance) to migrate to r, (nearest oir tane ٠٣٥٠ 


E ise the limits of Eq. 3 from r to xr. and ፲ to fy. We have 


፲ጊ 


Í 


previously noted that the field gradient charges very slowly beyond about 
imem 110m the anode wire (see Fig. 3). Therefore, if we select our rise 
time discriminator in such a manner that pulses which spend their full 
energy within the active volume enclosed by a Il cm radius.are rejected. 
iWiewcan see that the assumption regarding the linear dependence of 
mobility and E/P is not too serious a concept since the pulses which 
depart from this assumed linear dependence(at high E/P values) are 
automatically eliminated. As it turns out, the ''dead volume" of our” 
counter will extend to about 2.75 cm, as will be shown later. Pulses 
which originate in the "dead volume" and extend outwards will be treated 
in detail later, as will pulses which do the reverse. (Section II-B. 7) 

| electron production. From the ው we now have a means of 
calculating rise times for electrons that are located at various radial 
distances in the counter active volume. Wedo not, however, have any 


mrereation of where these electrons will be located--in fact there wil) be 


bons 


millions of electrons generated along the path length of a charged 
particle that is stopped in the gas medium. Let us now turn our attention| 


to the formation of these electrons in a gas medium. We have previously: 






! 


; TT 3 ፡ | 
discussed the various reactions that can occur with a "He neutron inter 
action, and have determined that only three of the five potential reactions | 

| 


are of interest in our counter. Neglecting the energy distributio me ii 





the three events for the moment, let us now determine the ranc ia 
the various particles in the gas medium. 


Friedlander and Kennedy“ 


present an empirical relationship for the 
determination of the ranges of charged particles in a gas mixture gue 


relationship, adapted for our specific usage, can be used to determine 


۱ 3 : 
the ranges of protons, tritons, deuterons, amd `He recoils as follows: 


Range formulas for the individual gases (mg / cm?) 


F TO (311: 
Rp in Kr (E) = Ro in air (E) [1.89 - 0.25 logjg (E)] - 0.36 1 
R» in He (E) = Rp iin air (E) [0.82 + 0.043 logig (E)] 


JHe Recoils 


| " V E E 
R35. IDEST UE) a ም inair (E) m" 89 - 0.25 1og1g (- GF 36) 
( 
73 ۱ ۱ E E! 
R35. in He (E)- 2 Rp تا‎ (E [0. 82 + 0.043 logig (=) 


“The energy distribution is discussed in Appendix C and in Section II- B. 6, 


mE 


Re i KEE) 


3 (Rp in air 5 1.89 - 0. 55 1610 (B)J-0. 34 
(6c 


) 
R4 in He (E) 


i r 2 
5 Rp in air 6 OT OO م108۰‎ (£z 
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ے جو0ہ 2 06 ۰(6 


Een ICE S) 


>) [1.89 - 0.25 logy ES ]- 0. T 
"(6d) 


à 2 {Rp inair | 
Ra in He (E) ta [0.82 + 0. 043 15510 172)|) 


ሬ ር l S IE 


The particle ranges can be converted into centimeters by use of the 


relationship: 


(mg/cm?)‏ مت 


መመር] Pie 4 (Z) Per 
14.7 14. 7 


EE CES 


where is the particle range in centimeters in the gas mixture, 
P = partial pressure of the individual gases (psi) 
P = gas density at STP 


ር particle range in the gas mixture (see below)‏ کت 


TT A ss (8) 
5:60). Rue (2) Rr ©) 


W = weight fraction, Wy, + Wx, =1 


Pije P He 
W ~ ا ات رر رت سو و و ہت‎ ር 
eo 88 PRI EE 
P 0 
Wye II KT = 1 ፦ Wie 


Pie PHe t PKr PKr 


where 
P = partial pressure of gas (psi) 
P = gas density at STP 
Pera = 0.16629 g/liter (7) 
፣ቲ ን... +) 


The above relationships enable us to conveniently compute the 
ranges ot the charged particle in the gas mixtures, and are amenable 
LO Tia CNIS computation. The ranges computed as above are plotted in 
Fig. 5 for comparison. Note the substantial difference in the length of 
the track for various particles of the same energy--particularly above 


about 1.0 MeV. It is this difference in track length that will enable 


us to remove the recoil spectrum from our spectrometer. 


4. Rise time discriminator values. The rise time 41٥٥۶ r 
for the pulses can now be determined since we ee a means of 
determining the pulse rise times and the length of track ot varr m: 
particles in the gas medium. I we recognize that it is only the radial 
component of the track length that is important in our determination of 
rise time descrimination, we can select tne most unfavorable recoil 
situation that may be encountered. Setting the discriminator Ds 
time appropriate to this most unfavorable ው we will remove m. 
consideration all pulses that can be attributed to recoils (He, .ںی‎ ۵1 
Xr) This worst case can be considered to be composed of two factors: 
(1) The maximum energy that can be imparted to a recoil nucleus, and 


(2) the longest rise time that this maximum energy recoil can cause, due 
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Fig. 5. Ranges of Various Charged Particles in Gas 
Mixture of 10 Atmos. Kr, 2 Atmos. 3He, and 
0. 5% CO, 


22s 


DE uu cowerucsl position in the counter. 

Appendix B (Case 1) summarizes the maximum energy that can be 
imparted to the various recoil nuciei. It can be readily seen that the 
3He nucleus can absorb up to a maximum of 3/4 the energy of the incider. 
neutron. The other possibilities for energy absorption are of considez =a) 
2-01۰۰۰۰2٣٣۷۷۰ Ihe Sle recoil is thus selected as the recoil 
٠۶3 و‎ be ranges ox the other recoil particles im the gas sS 
is likewise less for equal energies when compared with ^He. The worst 
CcOrmetriecal case that can occur is ii the track length of the recoil He 
nucleus is oriented in a directly radial po کے‎ with one end of the 
track at an infinitesimal distance from the wall T the cathodes This 
positioning is illustrated in Fig. 6. 

ise time discriminator setting can then be determine NEPE 
(1) calculating the maximum recoil energy that can be imparted 


Es. 0 9 EN (max) 


He 
(2) determining the path length in the gas medium, LE. (max) ee 
the equations presented in Part II-Bs 2. 

(3) determine maximum rise time by substituting into Eq. 5 the value 


b, the radius of the cathode, for re, and [5 -A(E3 (max)) | fos کی‎ I 
e 


. . . ግ ቆ . . . ! 
his rise time aiscriminator setting (t 


7 ) will then exclude all recoil: 


from the resultant spectra. 


mom 





3 | 
Fig. 6. He recoil track position to give maximum pulse 
የ ን 0 


5. Determination of radial positions of the nearest (r,) and furthest (rf) 





electrons from the anode for ere n,p)I and ?He(n, d)D reactions. 





Memcider the fact that to define the positioning of any event in space; itis 
essary to define six — of freedom.. Three of these degrees of 
freedom are necessary to define the positioning of the center of gravity, 
and three are required to explain the orientation of the event with respect 
to the center of gravity. For the purposes of discussion, letus ማችን ም ን 
that we can define the ''center of gravity'' as the position where the neutron 
strikes a ?He atom. Since the reaction occurs within the boundaries of 

the counting tube, we can define the positioning of the event (center of 
gravity) in cylindrical coordinates of radius, length, and angular dependence. 
The operation of the counter, however, is independent of the positioning 

of the event in terms of length (except for end effects, which will be 


discussed later) and completely independent of any angular relationship, 


። 


with respect to the anode wire. This means that the positioning of the 


5 O Cr V Ly can be defined in terms of a single parameter mi 





To define the orientation of the resultant ionization tracks produced | 
by 3He disintegration let us first turn our attention to the 3 Hen, p) x 
s nN FST በር 105 ሀሪ given energy Hand given recu ai proton 
energy En, we can determine the energy E of the triton particle and the | 


laboratory angle @, between the proton vector and the triton vector by 


use of the relationships developed in Appendix B. 


ጅጥ ዌ፡ ከ. ስፌ . ( 


wy 


and 


ረሳ.  ” ፐሜፐፖ E (1 


i 
243 Er Ep 





Note that& is always obtuse in the laboratory system of coordinates. 

Let us now define a coordinate system. The coordinate sys miu 
our calculations is selected as follows: 
(1) The event occurs at a perpendicular radial distance ro froni iE 
center line of the anode. The end point of rg defines the zero point 
S FO Oo. the Coordinate axis. 
(2) The z axis is selected to‘include the vector r and extensions thereof 
(3) The x axis extends parallei to the anode wire and passes through 


the point (0, 0, 0). 





xo d 


(4) The y axis passes through the point (0, 0, 0) defined as above, 


and is mutually perpendicular to the x and z axis. 


x 








anode wire 


Big. 7. Selection ot coordinate ys o 
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With this selection of coordinate axis, the event is angularly independent 
of the center anode . Calculations of the path m m of the 
ionizing tracks. will be made in terms of polar coordinates using these 
axes as reference lines. We must now determine the distances of the 
nearest and farthest electrons from the anode wire. For the purpose 
of illustration let us visualize the Ep and Km track lengths as solid lines. 
The orientation of these two vectors can have three degrees of freedom. 
Let us select the three rotational axes as shown in Fig. 8. 

Neglecting for the moment the ظط‎ vector, let us turn our attention 


to the Ep vector. 


We are able to define a position rz which we call the tip of the Ep 


TI 
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Rig. 8. Rotational axes. 


vector.. (We will refer to the Ep vector, although in reality we El 
speaking of a definite length of path of a proton of set energy in the. 
counter gas medium.) We can corde describe this POR r> with 
respect to our coordinate axis by the "n angle 05, the asmuthal angle 


22 and the length of the. vector 27 We can thus determine Io (Ep, 93. 22) 


Z 
1 E E 
92 > oP ፲ 2 p: 2, ?2) 
[ጋመ 
ر‎ x 
G. 
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p 
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We need an expression of ro in terms of rg since in determining migra 


tion times we are concerned with the radial distance from the central 


። 







Ep sin 8^ 
Ep cos 65 


Ep LI - sin? 6; cos2-g 2 


Ep sin @2sin(-¢,) 


Ep sin 02 65 (” ቀ.) 
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Plane containing x-axis ond vector E, 
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Fig. 10. Determination oí ፓ2(" 6 E, 92 92) 
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d diss 

counting wire to the point ro (Ep, 92. 05), i.e. we need an expression 

f ya e E 3 <= ደግመ ` +- ገ Fi - - 
Oi PIENE Ep» °. 95). To. determine this expression letus pass a,plane 
which contains the x axis and the vector Ej (8 plane can be defined by 
two intersecting lines). This plane will intersect the z axis at (0, 0, 0) 
and will make an angle 68 with the z axis in the minus y directii 
Rele rite rao TO we canere an Expression or ro(ro, Ep: 93› 05). 


By the law of cosines: 





2 2 
= r + E [ - sinf 9 ር68”(- 92) - 2 Yo 25:51] - sin 9 cos*(- 23) 


he 


| 
(7 cos (18€ 


but, cos (180- )=- cos 8 


TED CoS oO? 
Sd oI EIS IO cos > =. - 


ll 23H 
Ep 'l - sin? 92 cos” (- 92)] / 


substituting this value for cos B into the above equation, we find that 


2 £ 2 
95 = Ty T E - sin? 92 ones (- 92)] 202 ro Ecos 23 
or | (1 


በክ... 2 
r, (ro Ep 92, ሀ3)። 2 + BDL - sin” 92 cos TT Ep cos 05 


Let us now consider the possibility that since we have er placed a 
restriction on where our event can occur, that a skew line may Cs 
| 
and that the closest point to the anode wire may in fact lie somewhere 
along the length of vector Ep, rather than at the end point. Since Eg 
isa definite length in the gas medium for a given energy of the puo 


let us cali this length £ (Ep). An axiom of solid geometry is that the 


closest point of approach for skew lines is at a point where a line 


connecting the skew lines is perpendicular to each. This means that 

a minimum distance can be determined. We have an expression for the 
length r5. Therefore, if we take a partial derivative of ro with respect 
LO ES (note that this is actually X(Ep))and set it equal to zero, we can 
determine a distance Ao along the (Ep) vector whereby the skew lines 
are closest. We will cali the closest distance between the vector Ep 


SC the anode wire Ya. 


2 ; 2 4 y 
TAXES) l -= sin’ 92 cos (- 05) +2X(Ep) rp COS 32 ] 


|, / 
. peg - sin? 9 cos? (-ፆ#3) *Zrg COS 2 


Set the above expression equal to zero and solve for Zp; note that the 


e bracketed term can never equal zero, Forreal values of 145 we 


on tain 
; ro COS 92 
ዐ.ም ከ La 


where 0 a <J (Ep). 


The negative sign indicates that Üp cannot be positive unless 05 is 
ereater tnan 909. This is in accordance with the physical situation 
Sioned and the selection of the coordinate axes. Af 03 = x, p must 
equal rg and the Ep vector is perpendicular to the anode wire. ደፍ cannot 


exceed the physical length of Ep which has been assigned. 


RS a 


O C ተ ው this value of ¿p into Eg. ll for Ep determi 


Coat: 






2 T cos? E ጾ - sin? 92 cos? (- 62); - 2 ré cos? 92 






r. = /r : f. 24 . 7 
= 0 0 sin’ 0, cos? (- 0(7 |1 - Sin 92 cos* (- 22)] 
oa 
: cos? 92 
፣ ro - £ 


l- sinf 92 568” (- 22); 


where = <o, c 
2 
and 0 < < Z (Ep) 


Let us now turn our attention to the other portion of the ጋ (73. 531718 
Celts Let US Consider the ET secto 

In a manner similar to the determination of the end of the proton 
vector &p,we can position the end of the triton vector Erp. We shall 


call this position ri (Eq, 01, 01) 








wire 


Fig. ll. r,(Eq,9,, 01) | | 


In the case of the triton however, we cannot arbitrarily describe the 


position of the end oí the E» vector with respect to the anode wire, since 


there exists a unicue relationship for the positioning of the triton 
Eoctor for a given neutron and proton energy. (En and Epre nmen 
tively.) This relationship is characterized by the ancl 
the obtuse angle beiween the proton and triton. We must therefore 
Gonsider this constraint upon tne direction that the triton can travel. 


We can visualize this restriction by the following: 


2 
{ 
| 
| 
| 
| 





I 
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ANODE 





Fig. 12. Rotation of E. vector about Ep axis. 


If we were to rotate the Ep vector in Fig. 12 about its own axis 


while imagining that the Em vector is rigidly attached to the Ep vector 


4E 
at an angle ly trom the axis oí the E, vector, we can see that the tip 
of the E vector would trace a circle in a plane that is perpendicular 
ک5‎ > vector. The determination of the radial distance to the center 


wire ry(ro, ፻፡ጥ› 91. 01) must take this degree of freedom into account, 


Let us now examine Fig. 13. The interior angle of the cone described 


by the rotation of E. is (180 - œ ). Let us now pass a plane through tne 


“መ 
+ 


z axis and containing the E> vector. Since this plane contains the z axi 


۰.۱۱ 1610101601 be perpendicular to the plane containing the x and X 


| 
|] 
| 


axis. We will use this plane to relate 9, to 8», and 0) to 05. Letus 
consider the rotation of the E vector from a ''zero position!’ where the 
Em vector is at a maximum po sition away irom the center wire EE | 
POP position Par ንኒ ንሲችን።ን 7106117600117۴56 ٣۷ 2 1 NS Ep vector is 
then E in a counter clockwise manner (viewed from the tip of the 


Eo vector). The angle of rotation we cali $. By projecting the 








መጥ vector back to the plane containing Ep and .z we can deter the 


incremeñtal angle between 9, and (180 - 05) this incremental angle is 


n E sin cos» | dc sin Cl, COS y 
sin = ET 1 0 sin aaa 
|” S 8 sin 7) | | (1 - sin“ OQ; SEM 
therefore 0, is related to 52 by 
1 | sin Oy coS y 


91 
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180-05) & sin c MEAE E TES 
| 2) I 11 - sind Qr, sin y)! | 


OI 


| Sin cU COS. ^ 


(1 - sin“a L sin“ /፥ 


D 
ll 


"+8 E sud 


where CQ as "ያ. EN and En is a constant, y varie sm sS 


OtoZx. The sign ofthe incremental angle is positive for -%< y ረጃ an 
Z 2 


negative for L< y ERF ፡.”..ፕ፣ቭ ህጅ ፡ም ጅ ታው any selected. 
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oA 
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roiotion about ihe E» vector 
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Anode 
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Note: Vector Ey, shown exaggercted in size with 


respect to Cp for clarity of angular relationships. 
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Fig. 13. Determining Relationship Between 0> and 9 


1» and 
7 and 94 s 


If we likewise were to consider tne perpendicular projection of the 
LS Gof the کت‎ vector to the (ED, z) plane we can readily determi ٣ 
slnce this projection is perpendicular to the (Ep, Z) plane, it must of 
necessity be parallel to x - y piane. We may therefore compute tm 
incremental angle between @, and 9, and determine a relationship 


between them. Tne incremental angle is 


m 
E 


بے کا جو ری 
cos] | Eq (1 - sin^Q sin" >) / |‏ 


Therefore, the relationship between (, and 92 can be expressed as 


u ge 22) - cos” [1 - sin’ ay, sin? y] /? 
Or 
5] 


9; = 0^ T cos a - sin^Q y, sin“; yn 


where the negative sign of the square root is taken for 0% y < x and the 
positive square root is taken for x “jy <22 
C, is a constant determined by fixed Ep and En 
J» is arbitrary 
Having now established a relationship between 91 and 05, 01 and 0^ 
through the constantangle Ay and the rotational. angle y, let us now 
proceed to determine ry (rq, Ew, Io 91). 
Referring to Fig. l4, the same procedure can be followed as was 


used for the computation of 221" 6፡ Ep; 95, 02) in Fig. 10. We arz En 
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Fig. 14. Determination of TQ ጅጉ. 9); 


the relationship 
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Y. 7 f Tg m 14 - sin 91 cos (180 - #1) - 2 Yg E cos 01 ( 


Lhe possibility of a "skew" condition must likewise be conso. EEEE 
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in case of the Em vector. "۳ 4 ۵ (Ann) must be determined and 


substituted into Eq. 16. "The partial derivative of r, with respect to 
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KiE) must be taken and set equal to zero, 


Hos. 
ory በመ 2 m 2 
IE) = D z tAE) J = sin 9; cos (180 -$) - 2۵ Ern) cos ( 
T L 
ne "m 2 ; ] 
| 24024) f- sin 9, cos (180 - 04) - 2 ro cos 9, = 6 
i : 
Solving for £ we find that I 


pe TNCS 91 
E |l - sin? o ? (180 - 1) 
S 1 COS - 1) 


f, — 
TOT ED ሩዶጭ < AEN) 
substituting this value back into Eq. 16 and calling this closest point 


of approach 23› we find that 
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2 re E 9. =- S cos“ (180 - 91) | 2 Yg cos 94 

a => 4 = ። i. copre LO UE , 9 000 

3 እ ር 2 2 Enc f A 2 | 

ላ 3 > sin 01 cos (180 - Qj, 11 > 812 9, cos (180 - 01); 
D | : (18) 

ER h cos’ 91 
نر می سر ہر ےد کہ‎ a 
2 x | - sin 91 cos? (180 - 6) | 
burmumarizing the calculations for radial distances tron tice ህህ 
ome, given EN 5 "۵ we must caiculate four additional radial 


distances for each ms NUI CC Le UE 
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~+Q- Ep (9) 
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ro - 
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“4 L (1 = 613“ ፀ ። (- 0 ) ' 
A i EC 2 
(only necessary to compute r, where 2... ፓ 
and O«4p < Es) ۱ | 
Lg 
L D = ag 0 1 E سنہ ہہ‎ e NOE GI TIS = 
ላ 1 - 818“ 91 605 (180 - Ø; )| (18) 
and: Š : 
i ۲ Sin GT COS O . 20 ٦ | 
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The sign of the incremental angle is as follows: plus for - کر سے‎ 


a ۰ N 7 
mom St Toge?” 
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e ካምና 


0-05 T oci [1 - 5:6" ርኋ ፲ کے‎ 

The sign of the incremental angle is taken as negative for 0¢ 8 <7 
ana pesitive tor መበህርሪ/ ፍረ 

t nas thus been shown that all possible combinations of radial 
distances from the center wire can be calculated in terms of known 
En Ep, and rg in terms 02 05, 23. ፓ, anda. 


፣ 


Tí رم2‎ Ty, Tos رو‎ and r, are calculated for a given event, and the 
maximum and the minimum radial distances taken as the positions of 
the most distant electron, re, and the nearest electron x, to reach the 
central counting wire, we now have tne distances to substitute into the 
Gquation tor determining rise time. 
Deuteron radial distances 

For the reaction involving deuterons, there is an analogous m e m mE 

ና * . . > 

ship to the proceeding if we replace Ep with Ego, E> with ED? ande, 


with tne relationship: (see Appendix B) 
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P Fo | a - sin" 0, cos (- 2)] ro Ep2 cos 9, (21 
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cos” 9 
EN 0 1 ——v MÀ (22) 
4 |1 - sin 92 cos? (- 92)! 
(only necessary to compute Y, when -= 8 0ے‎ 
ء۰‎ p2) 
= 
; cos? 91 7 
::20፡ 0... - —— m A‘ CAERA) (e 
> : 1 - sin“ 9. cos (180 - 94). y Dl 
L l | 
I sinet, COS. 
and 0, = 180-0, + sin . (25) 
1 2 " E وہ‎ 2 
BE REED COT SUM > 
i : š ' SELAN s 
Ee sien of tne incremental angle is as follows: plus for سح‎ > JS : 
minus if a 10.2605 F -1 کے ےت‎ . 271! /፡ 26 
us 2 > and 0, = 02 cos | ٭‎ 388 1 (26) 


where the incremental angle is positive for ٭‎ %<2< and negative for 


OCT. 


6. Development of the probability function (Pp) for the replacement 





of pulses removed by rise 





We have now developed a means whereby the rise time of any valid @ 


pulse (i.e. not a recoil) can be determined for any possible energy distri 


tion and spatial configuration of the 3 He (n, p) T and the 2He(n, d)D reactidll 







۶ 
2 


If we were to vary Ep, ro, cos 9, Ø, and X over all possible ranges, 


O ° 
a tabulation oi rise times t, can be made for each incoming neutron ምይ 
Before tnis tabuiation can be accomplished however, three additional 
conditions must be considered. 

(1) The volume probability of an event occurring at a distance r is gre 
at points that are distant from the anode. In Zig. 15 below, anjemem s 


of the counter tube is shown. The voiume probability g(r) of a reaction 


Seen renecraeur Witnin a differential element dr is: 


where b is the radius of the cathode. 





Big. 15. Volume probability in the radial direction. 





ED will vary from Ep max tO Ep min for a set neutron energy. Any en 
within these limits is equally possible. Cos s selected as a variable J 
since the selection of Q would result in undue weighting of the rise times 

near the polar angles. A discussion of these factors is included in 
Appendix C. lt is assumed that the applied voltage is held constant in the 


Therefore, events occurring at various distancesr Ironi here nC S 
ihe صدوء‎ 6۶2< must be weighted in accordance with theiz volumeae o 
ability. 

(2) Thus faz,the finite size oí our tube has not been considered in the 
determination of rise times. It becomes obvious that some of the energy 
meecors will intersect the walls of the cylinder, and that the ionizing 
tracks wiil be stopped short of their full length. In this case a pulse of 
lower energy and rise time will result. An additional constraint must 


EUerore be placed upon the calculations of the various r's (1.6. Z r 


ው Ze 


T2, ry) in the preceding section, viz, if any of the r's exceeds b, then 
Be pulse must be considered as a "short pulse." Short pulses will be 
divided into two possibilities: 
(5). The radial component of the track lengths is still sufficient to cause 
| I 

the pulse rise time to be above the rise time discriminator level te This 
will cause a count to be registered in a lower energy channel. The rise time 
158 computed using b as r e. This type of pulse will be designated as a 
"wall effect pulse'' (W.E.P.) and will be totaled separately. 
(0) The radial component of the track lengths is insufficient to cause a 
meee tO be registered in any energy channel. The rise time will be 
computed however, using b as the maximum of the r's and catalogued in 
the appropriate rise time channel. 

for "bookkeeping" purposes itis necessary that the above two possi- 
Ties for "short pulses" be catalogued in the same energy channel 


appropriate to the n ©) considered for a pulse that goes full length 
፲ 















in tne track medium. They are not to be catalogued in the energy ehanned 





corresponding to their reduced length. Wail. efiects are separately E 
and if these degenerate pulses were removed twice, it would distort the 
ت2‎ 66 0371000 

Each time the rise time discriminator setting is changed the above 
two possibilities for "short puises'' must be recomputed. 
(3) To compare the rise time spectrum from the (n, p) and ٤۲ reac m 
it is necessary that the two reactions have a common reference conditions| 
This common reference condition can be met if itis stipulated that the 
number oi rise time situations considered for each incoming neut o mE mE 
be the same. The variables for the rise time spectrum (for a set neutron 
energy END fixed rise time discriminator setting, and fixed applied voltam 


are: three angles (cos 0, 0, and 7), distance ry and an energy distributic 


GES The angle variations are identical for each case considered, the 


variation of rg will be identical, but the = can vary from an E to 


P max 
| 


ማሜ. 


ያ : 7 ሽ 5 . 
Ep min’ the values of which are dependent on the energy of the incoming 


neutron. We must therefore consider in our tabulations that amjgdememeam 
number of energy increments (not identical energy increments) Os | 
between Ep max and Ep „jg regardless of the energy of the incoming 

neutron. Note that since only a ratio is eventually desired, this comori 
merely increases the numerator and denominator of the ratio constan 
E ctor. | 


All three oi the above conditions are likewise applicable to the deutezk 


reactions, as shown in the following discussion. 








(=ar 024 O = 
ዒ ባ ባ 
ዴ.-› 95፡01 £c pesi 2N/ £—. (pos ON 
I 
8 
q = q = + q = + 
Í 
t 
t- 
l: / 
I 
وی9‎ 4 ٦ SUIT} esta 
299870 oes? 
q 8111۸ ۹ 8111۸7 ۹ 
“ x ١١ 8و‎ £ d < 
کے‎ [n e Sos[n X 
jo | Hs 
(2osií Qc)N Qiequmw |^ * * S s 2)N iequunw | (29s DN 





(o + Nr) KÁS8xeus ye soU ƏSA 


——— —À 























በር 





— 


SnotieA 10; SosTnep- (du. Jo. vocum T 91i 





q=. 
7 q 
P; 
ጋ 
q 
ባጋ 
e 
7 
QUIT} S13 
ን) sosqnd 
UIA 29 VORBEI, 
8 ۰٤ TX TEM 
Jo Jo 20 (ito) 
iequanN | (THM)N | Toquiny 1 ና! 








a ee 








- سے‎ መመመ II ብ = _ a o 


n 


A tabulation of the number of occurrences of pulses having the 
same rise time for each (En + Q) channel can now be formed. This can 
be accomplished in the following manner {see Table I): For each 


(Ex 2:09 ) channel, i.e. the energy of the incoming neutron plus the 
m 3 








Qn.» Of the reaction, a reaction site r, is selected. The energy distribu 
ac 


tion of the particles produced by the reaction is varied in a numic mem 


-— 


۱ ا رج‎ E 
energy increments from Ep ,.4, to Ep Note that the total energy 


ees uae 
available is constant at (Ex, +Q) since we have not varied the energy of 
the incoming neutron. The tnree angular relations are also varied Over 
all possible configurations. For each energy increment, and for each 
angular configuration, the pulse rise time is calculated and tabulated. Tne 
number of pulses that have rise times at 1 Rsec, 2 usec, Cos 
summarized and multiplied by the weighting function n The ۹ ri 
site 1s moved a smali radial distance outward and the process is repeated; 
This tabulation continues until the reaction site is an infinitesimal distance 
trom the cathode. The volume weighted number of pulses of the same 


\ 


rise time is tnen totaled: 


[r-b Dm = o 
7 a 
s Ncounts'dr- N(WEP): Z + N(1 sec) = QN (5733586605. t. 
- OR 
C= r= ገ F = O 


| 


o 


H 


xls 


N(30 psec) = 1 
| 


1 
Or 


36 


wnere dr can be any constant and is added merely to remove the dimensional 


I 


ES | 
dependence of cm . For finite differences in r, the interval between the 


m" 





Er. and Yu COuid be selected, 


ihe apove tabulation is repeated for ali values of (Ex +Q) from 


e 


Ew 0, to S. 714.0 MeV. Note that because oí the previously imposed 
meme.t10ns, that the summation of all the counts for each (Ew tQ) tabula- 
pon are equal. 

Taoulating the 3He(n, p)T and the 2522 0:5, d)D reactions separately, 
Pyemwean Construct a matrix for each reaction which is a tabulation oz the 
number of events having identical rise times in each energy channel. 
such a matrix is shown in Table Il. 

An identical tabulation to Table Ii would be constructed for the 
(n, G) reaction. 

The (n, p) and (n, d) tabulations could be directly combined), ٦٣ 
reactions are not equally probable. However, if we multiply all tne 
entries of the (n, d) tabulation by a factor that is proportional to their 
relative probabilities, then we could directly add the corresponding 
elements in the two matrices to obtain a combined matrix of rise time 
vs. (Eu + Q). Such a factor is available to us in the form of the cross 
sections of the two reactions, oh, PUN) and on, a (Ew. (Note that these . 
are the cross sections of ፓ(ጅኢ) and not c (Ex + Q).) The combined matrix 
would appear as shown in Table III. 

Therefore the percentage of counts removed per energy channel is 


Paral to the sum of the counts to tne leit of the rise time cutoff (see 


Table Iii) divided by the total number of counts in that energy channel. 
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If the percentage of counts removed is designated as 'Pp,'' then the 






percentage of counts that were counted is (l- P4). Thererore, consicerAmi 


s ፈ . . INPS poe L ps ! | 
the raw spectrum that is received as output Írom the multichannel analy = 


and designating the counts from ary individual (Eyt Q) channel as Cx: we 


CN 


can say (1 E UE - CN- OT CAN E , Where Caw is the valid 


(1-25 
counts we would have received without rise time discrimination. 

- can be computed ior each (Eyn tQ) energy channel oí the multi- 
channel analyzer, and after correction we will obtain the corrected onc ٢ | 
spectrum with the rise time discrimination removed. 

7. Discussion of error introduced by assumption of linear depene sm 
of electron velocity and (E/P). 

it has previously been assumed (Section I1- 8.1) that the mobility of 
the electrons in the gas medium is a constant value and that the 90 ( 
of the electrons in their inwarG drift towards the anode varies directly 
as the (E/P) value, i.e. vzm» E/P. This relationship appears to be a 
reasonable approximation in the regions where tne field gradientes varies 
slowly. It can be seen from Fig. 3 that from tne cathode to sP sS 
distance of 1. 0 cm from the anode, this condition of a slowly van m 
(E/P) holds quite true. At distances closer than 1.0 cm from the anoaes 
the field gradient increases rapidly, becoming higher the closer the anode 
is approached, At regions very close to the anode, gas multiplication 


fakes place, and an ''eléctron avaianche'' is produced, 


The region of interest for electron migration time calculate 





nowever, is tne region between the radius wnere gas multiplication begins 


7) 


ረ 


(i.e. a iew mean tree paths of the anode--or as deseribed by Wilkinson, ` 
k distance írom the wire surface of the order of its diameter'') and the 
region where the field gradient begins its rapid assent (i.e. about 1.0 cm 
radius). 

First, it is known that the mobility of electrons is not a constant 
value over unlimited variations of E/P. At some value for eaeh eomponent 
fers, 2 Saturation eondition is reaehed whereby further increase in PE 
does not result in an inerease in eleetron adrift velocity~-and in fact, the 
veloeity may even decrease with inereasing I2/P. These saturation 
eifects are probably the result of eleetron collisions with gas moleeules 
EET re no longer purely elastic, and a portion of the eleetron's kinetie 
۴۰۰۰٠٠۰٦ imparted to the gas molecule in an inelastic eollision. Tne 
probability of an inelastic eollision apparently increases with increasing 
electron kinetie energy, but the inereased energy losses due to the inelastie 
collisions cause‘the drift and agitational velocities to remain at a constant 
(or nearly contant) value. In the ease of krypton, this terminal velocity is 
Beecned at an E/P value of 1.5 volts/cm - cm Hg 2 when the krypton has 


20 


1/29» CO, intermixed. English anda Hanna ` attribute this peak velocity 


ES result of the 'Ramsauer Effeet." On the othcr hand, no such 


۔ 


: ! : ۱ 24 
E inal velocity peak has been reported either in pure He, QNS 


SIS 


He- CO, MUTES., In addition, the veloeity vs. “@/P curves that 


t» 


are available from various experimenters do not go beyond an E/P vaiue 


-i € e ግ 20 
memere tian ©.0 volis/cm'- cm Hg for krypton; and 10.0 
E SIE 7 ۱ 
Col: Cm He for "He, whereas tne E/P velue varios 


- 1 


from 0.75 to over 190 volts/cm * cm ^ Hg in our counter, in the regiém 









under consideration. It is thus impossible (at present) to even estiméte 


the electron drift velocity in the inner 1.0 cm radius of the counter tube f 
| 


Quantitativo we can say that the electron velocity in this To E 
Y, y ር | 


will be considerably less than the value we nave assumed. This ohengzoE 
of smaller and slowiy changing velocity is due to the saturation erifecetsem . 
the Ar-CO, mixture--and probably an identical effect from the > de “CO, 
mixture that occurs at an E/P value higher 70 those reported. 

Let us now examine the eífects of this uncertainty of electron velog™ 
in this region in our counter tube. Equation 5 states that the ۷۶۰ 


a pulse is equal (for a specificd applied voltagc) to a constant times Bie 


2۰٠ت‎ 6 101 ۲6606008-268 Of two radii, i.e9 


ct 
1 
Q 
© 
1 
ርዕ 
e 
H 
1 
E 


We have likewise pointed out that for tracks that occur entirely within 
the region where E/P varies slowly, that we may consider the "rise tir 
as the time that it takes the furthest electron to migrate to ۶: ٤ | 
to that of the origina’ position of the nearest clectron. 

Let us now go back to the determination of our original setting of 
the rise time discriminator. The value of the rise time discrimination 


ہم 


! ጫ . 5 . ° 
t. was determined by solving Eq. 5a with the values r, - b and 


^d 


አክ 5 à 3 ae ግ ግ x 
riz - Z... (= Ex). Once this value has been determined, we can 
TG Js 


interpret Eq. 5a as the- area between two concentric eiL C E 


í 





4 1 
T 


ig. 16. interpretation of rise time values as a constant area. 


The total area between the circles remains a constant, i.e. 


2 Z 1 D 
2 C SUAE 

V. £n ^c 

| D 

ان 

wnen r, = O0, this expression reduces to: 
a Zn 
reina) 7 | (55) 


Where T t(min) is the outer radius of tne "dead volume, ' and any tracks 

۹ u . . . A . + ግ 

Peewee completely. within the cylinder bounded by Tf(min): ፣" ባሻ. 

መክ ከርሮ by rise time discrimination. From the above discussion we 

ከ )25 that Tf(min) is not as large as calcuiated in Eq. 5b, due to the 
۳٣۰۰۰۰٠٦٦٦٠٦۰۶ ٠ط‎ electron migration velocity in the region x. ۷ءء تر‎ 
since the velocity tends to be considerably less than the linear relation- 


snip originally assumed in a nigh E/P region. For ai14.0 MeV neutron, 


. ግ " = . UR. . . q 
whicn procuces a 10.5 MeV maximum "re recoil, we find that 





i 


+ = 6.1] usec and rz 


: fGmin) ^ 2. 7 cm (for Vo = 5500 V). Since rf(rain) 


(actual) is less than the value calculated for a fixed t the ''dead volume" 
of our counter is less than calculated. The question of 'how much less" 
is one that cannot at present be analytically solved. It has been 
oreviously stated that it was considered that the linear dependence of 
velocity with tne (E/P) factor was valid at regions r>1.0 cm. If the 
Yf(min) (actual) is not less than 1.0 cm, then our assumptions arcas a 
valid for those events whose tracks lie entirely in the active volume 
beyond 1.0 cm zadius. 


We have assumed that the probability o: a reaction occurring in any 


unit volume of the gas medium is identical. The volume enclosed by the 





n 
1.0 cm radius comprises only 4% of the total sensitive volume. The track 
orientations can likewise be completely at random, and so the probability 
of a track having any portion of its length in the volume enclosed by a 

1.0 cm radius is dependent upon the disintegration particle energies ummm. 
reaction site, and the orientation of the reaction particles. While the 
probability of passing into this "inner wall!" region can be computed in a 
manner similar to the probability of a particle striking the 1+07 

(see Section IT-C), there is no way to relate this probability ٦ 
time spectrum, since electron velocities are unknown in the regio mm 

also every accepted count will be registered as a "full energy count." 

In addition, even if the tracks of the charged particles do enter into the 


region of the "dead volume, '' their total radial orientation may be such 





ፎሮ 


that the pulse vise time is greater than t ana tne counts will be 

it can thus be concluded that: 
(a) Pulses that are the result of tracks iying completely inside the "dead 
wOlurme' will be excluded due to rise time discrimination, regardless of 
energy and orientation. (4976 of the active volume.) 
(b) Tracks that lie completely outside the "dead volume" may, Or may 
DDE counted due to rise time discrimination. Tnose pulses apove a 
" ን rum energy that were سے‎ aque to track orientation can be 
Mathematically replaced. | 
(c) Those tracks that have some portion of tne length within the ''dead 
volume'' may, or may not be counted. Tae accualo rise tired S oe 
purses will be larger than the calculated rise time, but the rise time 
discriminator setting will be the same as in (a) and (b) above. This means 
that more counts of each energy would be registered than would be indicated 
by Our theory (i.e. more counts because more pulses would have longer 
Bese times and thus more would escape rejection due to rise time discrimina- 
HOn.) 

The magnitude of the error introducec by (c) above is unknown, andits 
effects must await the comparison of the final spectrum as developed from 


Sm Counter with a spectrum obtained by other means from the same source. 


E Dotanticipated that this error will be large. 


persons who are interested in development of these functions, but in 


Eros 


i DEVELOPMENT OF A PROBABILITY FUNCTION FOR 
THE UNFOLDING OF THE WALL EFFECT 






Te development ol a probability function to remove ۶۳ effects wi 
primarily developed by Brown, end then extended by Wang. The overall x 
development is quite long, being in general the substance of two theses, 
and only certain aspects a be presented herein along with 8 ٤٤۶۶٦ 
to Wangs "Wall.Effects Function". For a complete derivation of the entire 
Series Ox provabilivy ር E are referred to Ref. 9 and 16. 


The notation of Brown and Wang will be continued in this paper to aid 


general the..: functions will only be presented and identified in Ci P 


final Tom.: 


L. Wall Effects 


Tne sensitive volume of the proportional counter tube was previously Jf 


I 


identified in Fig. l. Within this sensitive volume of lengtb DEM | 


radius b, let us select a reaction site at r, . and assume that the 


Í 


| 


track length of charged particles is E . For the analysis, them uem 
gas filled volume (not just the sensitive volume) may be considered to be 
divided into four regions which will be identified in our ss 
discussion. 

à. Region A 


Region A is the core of the sensitive volume defined by raáius Of 


b -Ê anā length L - 2 Cas show in Fig. i 
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T's ለግ ነ ግኀኣ ገ 7 Dmm ለ 
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It is clear that only tracks of full size can be found in this 
region. Therefore, the probability PCP) per unit track length that 
a track will have a 


length between @ and Ọ + dQ in this region is: 


፻, (ዎ)ፅ፻ ። 6 (ሮኖ - ይ )ፈፆቦ 


For a single reaction, happening in‘region A + the probability is 


OO 


P, (PaF = | ድፎ -ፊ)፻ » 1 Cr) 
O 0 


Oe Region B 


 . Region B is the hollow cylindrical volume with an inner radius of 


b D an outer radius of b, anda length L -= M in Fig. 18 


CC 





Figure 18  - Region B 


For region B, Wang classifies the track length into two categories; 


Case l 


For a single reaction In this region, those tracks not intersecting 
ምሪ 
the wall will have a Tenetn ae | the full sized track. The probavility, 


pue), per unit track length that a track will have a length between 


© and 0 r d OÓ) in this region is. 


P CAL -=r EEE -L (t) 


where the factor Pul (P) represents the fractional arca of a sphere with 


eae 


radius 4 centered at (0,0,0) that falls inside the sensitive volume. 
Referring to Fig. 19, it is seen that the fractional surface area 
5.۱1676 136106 Lhe cylinder is 


un 12 [2 sinoaoag 


7: (ቦ ) = —— —— = l- 5 / sin 0 à 9| 
lx ° 


ግ 


wnere S = sphericel surface area ouvside of the cylinder. 


Brown makes an approximation that the transition -from 9 to 9, along 


by the sin“ Ø, and he.concludes that: 
ሯ 
A 


counter wall 





counting Wire 


Figure 19 Intersection of a Sphere and a Cylinder 
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Case 2. 





Tse tracks that do interest the detector wall will lose a portion 
of their length. The probability of ‘hitting the wall is then the spherical 
surface outside the cylinder i.e. 


25-5952 
1 |P ° - O i b -2\-3 
Pyo(t) = = 5 E z (2 25. e 


Then the variation over 3۳9.٣ ee f 
B2 eoa? 


Ads 
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c Region C 


Regions C are the portions of the sensitive volume of distance Q 


irom the ends of the sensitive volume, as in Fig. 20 





Figure 20 - Region C 


The track can be classified into two different cases. 


Case 1 
For a simple reaction in this region, the probability of having a 


n ek Length is 
20; (ር) $e = Foy (e) & (£-e) ae 
E ሮ |. مم مه ە مده‎ = 1 - 2 Join ህጻ 


ብ ም ሎ፦ት፦። ሙ፦ 





Figure zl - Region C 


Inversection of Sphere 
with Active Volume 





The probability of having a curtailed track length is 


DESC (1 sin a da) xe 


27006-٦6667 بی‎ = Za, sin 'O-dg = 22 cules therefore 


P. (Q)ae- 22 9 if z < تم‎ 


G2 ( 
= 0 ir z > O 
C. Region D 
Region D &re the portions of the gas filled volume beyond the 


Sensrcive volumes, as iN EIEN . 





Figure 22 = Region D 


Nuclear reactions occurring at sites outside the sensitive ۶٦٦ 
can also cause energy to be deposited inside the sensitive volume of the 
counter. 

If the procability of depositing energy from @ track length inside 


the sensitive volume per unit track Length of P is PPM then: 








። 


Py E Z sina da x2 


DUT COS Q 


l 


pnereiore sin a dO = = Uv) a for Lap and 


z ae E T 


^D (:)67 = [መር for ( ee z)2070 





meee Sens iti ye volume 
Figure -25-Region D; Track Lengths Inside Sensitive Volume 
. & Summary of the General Wall Analysis 


For an arbitrery track length £4 the gas-filled volume can be 


divided into four regions, as in Fig. 2h 


> 
S 
ATI 

مم 








Figure 2l: Summary of Counter Regions 


mec ۶000071716166 per unit track length that a Crack u i) n ee 


length between e and (e + de ) jin each region are listed in Table IV. 


ee ee ee - 
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f. The probability, P(E-5E.), per unit energy that reaction products of 
دج‎ 


inside the sensitive volume due to the 





Enerov E will deposit energy 
wall effect 

According to the foregoing analysis, now, if we assume that the pulse 

height is proportional to the energy deposited inside the sensitive volume, 


then it is reasonable to assume, for the first-order approximation, that 


= K, 

E; = Kp, 
By, = Kb, 
BE, = KL, 


where K is some proportionality constant that will be cancelled out shortly. 


\ 


Therefore 
P(E -»E;) = K'p(4 ^ p), 


where K' is another proportionality cónstant. For the whole counter 


Des 


Mp 
| Pp2(p, z)g2(z)dz. 
0 


From Subsec. II-D. lb. 


b 
p 
pot ج‎ of Pp2(p,r)gi(r)dr * U Pc2 p. z)g2(z)dz 
E 
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EE b ای وج کے کہ‎ B Z A 
Pnz: Eie E z ez S 9132 ” ይ. 


ሺ ክክ መሮ volume fractions are, respectively, 


2rd 
«1 (፡65 2፡1 ) پ9‎ ۶٦ 


dz 
g5(z) z, g3(z) 2 L 
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Therefore: 
b 2 5 
p (š T) = z / /ኔ TRE 1 + 28 . ጋጀ S cs 3 
/ W r TEM. :2 
D EC 
° 2-ዕ 
፣ 66 7٦7 ر2‎ 
f "nc — 
um (4-e + 
0 ዕ 
0 L ጊ i 1. 
Ja = لے‎ x ድ 056 7 
p» (£-e) ግ ከከ ጋሙ ከ 


(2ፐ) 


Where ES is the energy width of the particle in question, amg E. is 


the energy length. Note that this is the probability or fraction of 
iU 2 ከ. ር) ባተ ምታ” B, due to the wall effect. Surprisingly ٤۶٣ 
probability turns out to be constant independent of the energy of the 


| 


| 


incoming neutron.* 


2.  Probab&dlity of a Track Not Hitting the Container Envelope 


Brown and Wang derive expressions for the average laboratory energy 
of the various particles in a plane perpendicular to the counting wire. 
The expressions cannot be analytically resolved, so numerical Ios 


is performed on these expressions and the resultant energics are cics 


T 


approximated by use of the following formula: 





ر 
sis‏ 


E 
E, - 0.505 E. 4 0.450 


* WangS more complicated probubility function is aue to the improper 
selection of limits of integration ror region C and the approximation 
utilized.in determination of Pu CP) de. 
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COE 


td 
HI 


0.250 E 


y * 9.150 | | SS 


t-i 
i 
II 
A 
-i. 
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0.52 ፳. - 1.27 
ፈዛ 


tj 


ፌ 
pere the Cesignation E indicates. sre average laboratory energy 


in & plene perpendicular to ihe anodc resu ting from a neutron of energy 


E.. All units are in MeV. 


N 
From these expressions of particle energy, the 730617 ۹٦ 


Simultaneously emitted proton and triton, or a pair of 8606 or  + + ٤ 
not hit the walls of the containing vessel are derived. Using one 


notation of Wang, 
f. (E) = probability that a proton-triton will not hit the wall 


T(E) = probability that the deuterons will not hit the wall 


2 
ኋ (X - (۶2 
£, (E) - fE) : f. (E) m (29) "S 
22 
where f» (E) = the probability that a proton track will not hit 
the wall* 


flm = the probability that a triton track will not hit the wall 


2 
(1 - =a factor added to account for the forward track 


component of tne triton. 
JE 


CESTA, 1⁄2 where E 0.55 (Mev) 


: ( By i ሺ sin! 2 (29) (ተ) 


Fb 
~ 
td 
p 
Ld 
E 
— 
DR * 
H 


A 
* Both Brown and Wang calculate two different values for f,_(E) and 
equate them for energy boundary conditions 9ت‎ ۲06 ی٣‎ 20 r 
case, the energy boundary conditiuus (vecause of a larger size tube 
and greater filling pressure) are such that we do not ever reach the 


"conditions that Brown and Weng encountered. (i. e. the E..., for our 


| i T nia : یسوم‎ 
tube is at 17.0 MeV which is above the designed maximum ói" '1h.O MeV) 


ር 


/ um. | ds 
where ۰ = lencch of proton track in the gas mixture Zor a sr 066 


بد 
ንን ከ ከ P,‏ 


aad $5 = Jonc sl OC triton track in the gas mixcure for 3 ٦ 
of cnergy D : ۱ 
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£.(E) >: h _ i gj" (1 - 1 J | ቺ 


eee 7۶٣ 2 
"D (1 + =) 2 


7 E 
ም N To -— : 
and xy - ( جج جح‎ ) p 54 = 3.24 (MeV) 
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oa = length of a deuteron track in the gas mixture for a deuteron 


fon ener yE 


D 


The above nomenclature and designation will de. used in subsequence 


N 


sections. 





D. Combination of the Two Probability Factors to Reconstruct the 


O 


"e = Ts ር 2a 
2776110861 7660 7 ٦ 
RN E te II mw 


The original purpose of our neutron spectrometer was to accurately 
۳١۰۰٣73 a neutron spectrum in the energy range 77631" 6 8 6 61 ء١۷١5‎ ب٦‎ 
to 14.0 MeV. The spectrometer operation was such however, that in order 


to remove the more nunerous recoil pulses which were smeared across the 


Te 
و‎ 


energy spectrum from (2 N)... to zero, advantage was taken of the 


max 

differences in pulse rise times between valid pulses and recoil events. 

The removal of the recoil pulses, however, was not without its disadvantages. 

4 was discovered that to obtain a neutron spectra without recoils, it is 
necessary to forgo all information at the lower energies since the track 
"5 Of Lower energy pulses were shorter than the maximum energy recom 
unn These Lower energy pulses were thus removed with the recoils. m 
adéition, the orientation of the tracks of a valid event may be sucn that 
the radial components are Less than the minimum required to produce a 
ከ ከ፡።ን Liat exceeds the rise time minimum cut oft. 

The loss of the lower energies of the neutron spectrum is unfortunate, 
and & means to overcome this deficiency has yet to be devised. A means 
to effect the replacement of the higher energy pulses that were removed 
due to their specific orientations within the counting tube has becn 
developed and presented. earlier in this report. The probability function 
E) can be directly &pplied to the raw spectrum fróm the multichannel 
analyzer to replace "rise time removed" pulses in this situation. 
Tne other portion of the problem in the determination of a reconstructed 

neutron spectrum has to do with the removal of the wall and end effects. 
As mentioned earlier, this problem was essentially Solved by Brown ana 
extended by Wang. Modification to Wang's probability function have been 


made, and the revised probability function E has been previously 





presented. 
| What remains to be accomplished is to present a computer program 
outline of the means to unfold the various effects and to arrive av a 
Pemplete neutron spectrum. 

For the purposes of this section and to add clarity to the computer 
program, it is considered that the following subroutines have already 


been accomplishec; 


D 


(a2) Assume that the raw spectrum has been corrected for rise time 
discrimination. The raw counts in each (X. + Q) channel have been ' 

L 
ES ECE BT une prooability function (3-2, ond the revised channel came 


are now designated C... 
8 e RN 


(b) The counter efficiency has been calculated, and an enersy dependent 
value of the efficiency can be shown in the various reaction cross- 
sections. The eificiency of the counter was developed by Wang to ٣ 


function of (a) the nuclear cross sections, (b) the counter dimensions, 


and (c) the gas pressures. 


(R) = N. G, : | Ag 
Tl np a^ 7 PS. nyp(EQD fQ(EQ) 

+01 6 z E 

N, = number of “He atoms/cm” of sensitive volume at filling 

He 


- pressure 


60 te (n,p)T cross section (see Fig. 2) 


"n.p y’ P 


UJ 
| 


- Effective thickness of tho Getector for 4x geometry 


Total sensitive volume 
effective area 


where sensitive volume of sphere s sensitive volume of cy 


zxE = x b° L 


or 


(2 p +. 


/ * "NE 
thus efiective &rea - x n = e 


9. 1) = b 1, m T 


Le EE 12/3 
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| 
a رھ‎ = Probability the proton or triton ۹ ٣٦ hit the wall. 


ኃ 


A E expression can be written ior the efiiciency 2205 ٦ 


(n,d)D reaction 


Hn n,a Ey) 2) Id 


| ٦ب‏ ی۷۶۷ 
x pcc ۱ |‏ 





aes 


where the notation is similar to the above 


It can be seen that the efficiency can be reflected in a new cros 


section notation such that: 


where کے‎ a is the cross section thut has been corrected for the 
22 í 
eificiency of the counter. 


anc 


where o ns is the cross section of the Ass (n,8)D reaction that has 
| ; | 


[meen counted for the physical efficiency of the counter. 

Assume now that the higher energy rise time discrimination removal 
counts have been replaced, and the cross sections have been corrected 
for the counter efficiency as indicated above. A computer program to 
unfold the wali effects and to calculate the resultant true neutron 


spectrum could be 5 ۲0116151 


l) Let E., =” B, - | 
2 iE j Shp 
where کے‎ = Actual neutron energy 
E, = ‘Registered neutron energy 
(2) Let As = Number of neutrons of energy Es: 
ue E. 
3 j 
Iu Let Cay = Number of counts registered in top energy channel of 
the multichannel analyzer (after correction for rise- 
discrimination) 
E 2o) E 
(5) Solve for Ayr 
A URN 


= E x< E KA 
(6) ات‎ Acye1)! 9p Eg-1)! f Xn Ant pte? 


et, (Ey) | ° Pp (By "E. o) 





where P. (E -SE ) is the probability that an (n,p) event from 
ን N (N-1) 1 በ6) 

1 regi j annel EB E — E. ) = ج‎ 
channel E. will Es registered in channel (x-1) P( N E, ) بعد‎ ٦ 
وت‎ Solve TOY S 
(8) Repeat this procedure until we reach channel Aty) 

where: 

i=j-Ll 
= = à e ' د‎ (> — 
“(5-4)፤ 7 =05-4) ” -ር5 کے‎ EG) P(r) | اجچھاہ‎ 
2 E کت‎ NS *(፳-4)' 
(9) Let ت‎ ን ES E and the deuterons begin to appear. 
P l Xx 
(0) %0-4-3) ” *(8-3-3)፤ %.ጩኞርኡ4-3)፤” fiie)! * (eas 
| i=0 


Gi) E (Eeri) Pp Eaa Fi-a) + Ane ار یھو‎ 


f 5 Uu. E -1)| The last term in this expression indicates the 


deuteron effects from E, anā f (E ) = the probability that the 


N' Cee!) 


deuterons of energy E Will not hit the wall. 


(N-j-1) 


Dr Solve for A 
( ) (X -j- PED 


L= j ¬ 
(8-3-3)፤ ” %ክ.4-3) ” ነ *ቤ-4)፤ %ጩሜቤ-4)፤” 1-5 66.4)" 


=0 
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r. 


Ko چب‎ DU A m a n) ٦٦ 
divided by ፎ (m, . x 
۱ | E رہ‎ m ىك‎ 


ኣ 





(lz) Then in general: 
| 4 =.44+ 
ሳ፳-4-3-5)፤ ከ ኤኔ ከመ መሙ کر‎ | E Tu SC Tas 


سے 


E 8:64) ° Pp NY Png-1-k)^ 7 PGexk)s ٹک‎ 


j= kaem 

fs (Eg) ነ. Arai aa F(a) 6 Fars) 
ic x: 

ا یرت کے تا Pa‏ 


divided by p E- و‎ Do ii ظا‎ 1x) 

- Ea . "r a S 
where P, وت و‎ is the probability of a deuteron 
having scattering from a higher energy channel into channel መን 
ana f< ات‎ = the probability that the deuterons having energy 


5(8-4-3-:) Will intersect the wall = [i-t .لی‎ 


(15) This procedure is followed until the spectrum is unfolded to the 


lowest energy channel that has escaped rise time discrimination. 
(1i) The unfolded spectrum 


Ay Sed Atg-2)!7 را‎ E E کے‎ 


is then the neutron spectrum corrected for (a) rise time discrimination, 


(b) wall effects, and (c) counter efficiency. 


$ szi 
- DS 


` 
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+ 
, 
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removal of the rise time discrimination effects, the wall and end effects, 


ےد 


E. Random Probability Inversion Matri- . 


The foregoing sections. have considered in detail the computation’s 
necessary to change the raw energy spectra to a neutron spectra by the 
Euge rrecHuonssorethesspectrornete ۰٠٢۰٦٢ ٠ 

Unfortunately, as can also be seen from the foregoing sectione Emin 
number of computations necessary to remove these various effects and 


corrections is exceedingly large. "There are four position variable mM | 





ላ 


independent energy variables (one reaction particle variable, one. incoming 
neutron variable), risetime discriminator setting, and an applied voltage 
condition--a total of 8 variables. The number of variables and constraints 


on individual computations make their solution by analytic means virually | 
^ 7 | | 


an impossibility, and even when considered for high speed computing 
machines, the number of computations involved would require several ] 
hours of machine time. It was therefore decided that in order to compute 
an "inver sion matrix" which could be applied to the ''raw spectra" output 
of the multichannel analyzer recourse must be made to a "random prob- 
ability method béing that it will require only a comparatively small 
amount of electronic computing No un. and in addition will permit | 
computation of all corrections to the raw spectra simultaneously, as well | 
as elimination of several assumptions that previously were necessary to 
apply to the more general case. 


The random probability analysis for our counter has been developed 





ም 


by means of a FORTRAN PROGRAM having the code name "CATMAN. " 
r outline of the random probability analysis is presented below: 

ELA neutron of à specific ን En (selected as the ۰07٠٥٠٣ 
energy that our counter is expected to detect) comes 7 counting 
tube. The direction of the incoming neutron is selected at random, and 
since the detector is to have a 47 re sponse, this direction s spa 
ss pletely arbitrary. 

(2) The neutron passes into the volume of the tube (not necessary to limit 
this to the "active volume'' between the field tube S) and transits a distance 
into the gas medium. Since each unit volume is presumed to have an equal 

pusbability of à reaction occurring in that volume, the position of the event 
I elected at random by a ‘random number generator’ Which deter mame 
both تحت‎ distance along the anode wire, and radial position rg. | 
(3) An ہے‎ occum at 4A TAdIUS 7 the counter center line and at 

: | à 

some length £ along the counter length. The event can be either a Je e 
reaction, ~ a 2276 (9, 9): reaction. The probability of the type of reaction 
is based on the ratio of the cross sections at the energy of the incoming 
neutron to the sum of their cross sections, i.e. the probability of a 
3He(n, p)T reaction is: 


On P (EN) | 
9n. p VEN! am con d Ew! 


and the probability of >a. He (n, d)D reaction is: 


T HUC A Eg í 


= . 
2 /ሥ 


7 
۰ 
میود ر ے 


(4) The reaction particles are assumed to have isotropic distribution in 


the center of mass system of coordinates, and the direction of one of the 


-፦ 





réaction particles is selected at random (three dimensional selection 

based on 6j dcos T and equal but opposite momentum vectors). Them 
ር | 

energy of the particles in the CM. system is divided among the particles | 


as follows: 


EN 
Ep CILE (ZEN Qn, p? 
ሩ 4 3 . 
PE m He(n, p) reaction 
“Teme. SD 
E = E . 5 + Q ) x. (n, d)D reaction 
Bc ሪመርገን ጋሐ N ma ሽጋ res ከበ] 


(5) The reaction 3s converted into the Laboratory system of coordinates, 
and the energy of the various particles is computed in the Labor ë 
system. Note that any energy between a maximum and a minimum 
Laboratory energy is equally probable and can occur for the proton and 
ea eon particle 58. Likewise, if ZEND i a| (threshold energy), the 
two deuteron particles will vary between a maximum and minimum 
eons ory energy with equal probabilities. (From Appendix Gp 
Laboratory E. are determined by the random selection of Bem and 
5 O oa: 
(6) The length of track of the various Laboratory particle energies are 
comiputed by the relationships previously pre sented. 

| 


(7) The track lengths and orientations are then compared to the reaction 


site to determine if any portion of the tracks has intercepted the walls 
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of the tube, and/or has a portion of its length in the gas volume determined 
by the perpendicular radial projections of the interior ends of the field 
tubes and the end walls (end effects). 

(8) If the track of one of the particles intersects the wall or the e 
portion of the counter, the portion of the energy represented by the track 
that is physically within the active volume is calculated. Note that this 
٣۰٦٠٦٠٦۰-۶ 1s not directly proportional to the track length since the emer 
deposition is not linear (see Fig. 5). The energy deposited in the active 
volume is equal to the total energy of the particle minus the energy 
represented by the length of track that passes outside the active volume 
(if the same gas mixture were present outside the container). This 
conclusion can be drawn from the concept that after a charged particle 


መሙ” 


has lost a certain amount of energy, it still has the same length of path 
to travel as an identical particle that its life with an energy equal 
K reduced energy of the particle in question. 

(9) The various کت‎ values (radial components of track lengths) 
EN  iculated to determine rise times, s condition that an appro- 
priate constraint be placed on those tracks that intersect the walls of the 
e Gr pass into the end regions. 


(10) The rise time of the pulse is then determined, based on the time that 


it takes for the nearest (r,,) and most distant electron (rg) to reach the. 


i 


anode. 
IP 726۱-3 9 
፲ 2ጩጪጩዥኛ ٢ =n T EL "f n 


0 


(11) A weighting frunction is applied to the result to account for the 


greater volume probability at points more distant from the anode. 


ከ 


=— dr 
D 


m rss 
where r is the reaction site indeed in (3) above. 
(12) The weighted pulse is catalogued by energy (Ex tQ), rise time, 
and reduced energy (if applicable). | 
(13) Steps 1 through 12 are repeated several thousand times for the 6 
incoming neutron energy. The probability of the various particles 
Tino the voll IS computed. E the probability of the proton-triton 
tracks-not hitting the wall, andi.(E), the probability of the deuteron 
combinations not hitting the wall, are calculated after "weighting" for 
volume ہت‎ cc.isiderations. The probability of causing ampri 
a lower energy channel P(E —~E,) is likewise calculated from the track 


length-energy relationship described in step 8, after the volume) prem 





ability function has been applied. All three of these probability functions 
are combined into a new probability function P (E — F) forn spectrum 
unfolding procedures. 

(14) The energy oí the incoming neutron is reduced by one energy unit, 
and steps 1 through 13 repeated. 

02 The energy 5 the incoming neutron is reduced in successive steps 
until the owe st.value oí interest is obtained, or until the length of track 
c r 1376604 particles are less than the length required to 7ة‎ 


time discrimination. 


اوت 


(16) The weighted number of counts in their respective energy channels 
are totaled and accounted for, as explained in the preceding section. A 
matrix arrangement can then be compiled and corrected for the efficiency 


of the counter. The matrix arrangement could appear as: 











m K 7 
| Cy | » 1. 72 85 . 5 5 5 5 ° ° ° ain A1 
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| 
LE 
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| | | 
| 
E | ay ee AL 
Counts regis- aij = fraction of neutrons of energy Total neutrons 
tered in each Ji ጋ”ር counted at cme ug a sent into each 
2۰٦٠ incre- energy range. 


mehi for random 
probability 


analysis. 
(17) The problem now resolves itself into determining an inversion matrix 
whereby we can calculate the original neutron spectra from the counts 


ፌ 


registered in each energy channel of the multichannel analyzer, i.e. to 


‘compute an inversion matrix such that: 


ሥ 
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en. — 7 == a 

Ay by 1 bi ም) by 3 e ቄዊ 9 9 ቄ ظط ¢ © 9 ھی‎ Ci | 

A, b.) 522 523... 554 C, 
us | £ 
| = | x 
| | 
PA. | | - 
| | | CN-1) 

| ! 
| Aw | | Dd p ظط‎ e 65 9 5 6 6 6 ቆ o o ፡ | | Cx | 
= ھی سے ارے‎ L 
Neutron | Inversion Raw Counts 
Spectra | Matrix from 
| Multichannel 
° Analyzer 


(18) Note. that in stéps 16 and 17 we have referred to a matrix arrange- 
ment and an inversion matrix for our illustration. The actual matrices, 
if they were to be constructed, would be an array of 500x500 or more, 
and would be difficult, if not impo ssible, to invert. Instead of actually 
constructing these matrices, the computer would solve this problem by 


a subtraction-division process similar to that explained in the preceding | 


section II (D). 
(19) The complete random probability analysis in FORTRAN language 


(with appropriate comment cards) is included in Appendix E. 


ግሜ 











on 
IMI. CONCLUSIONS AND PRESENT STATUS OR PIOTR 


At the present time, the MES neutron spectrometer describe San ia 
as under development, has not been tested. The electronic d for 
the rise time discrimination however has been constructed ም operated, 
using the small site tube used by Wang. This tube is filled with 10 atmos- 
pheres of 7 (no heavy stopping gas), has only a 15/16 in. inside | 


diameter, and an energy width E, of only 0.5 MeV and an energy length 


16 


Er Eu 5. VeV "for protons. * The intrinsic tube efficiency (due to the 


small size) Ee ever iE extremely low, even at these moderate energies, 
and it was not considered that the expense of developing an inver sion 
matrix was warranted. | 

The "inversion matrix" for "ere larger tube that is on order ironi 
ESS Nuclear I has essentially been completed, and the 
"random probability program" (CATMAN) which is used to generate the 
iEn son matrix" has been completed. 

Future work in the development of this He SPECtrOmeter ۶٢ 
in some means of electronically comparing the rise time of pulses to the 
pulse magnitude (energy) to-determine if the pulse should be "gated'' into 
the multichannel analyzer. This comparison feature would be tm ٦٠٦ 


! : : a . . i : : | : 
Eur value (t ) for the rise time discriminator circuitry, and would 
T 


enable the registering of those valid low energy pulses which are now 


"These figures have been recalculated in this report, since Wang made 
an error in the density of ?He on pages 87 et seq. Oi Kerr p ss 


"o= 


automatically eliminated because their short track lengths cannot over- 
come the rise time discriminator setting. The rise time setting would 
then be some function of energy, i.e. t Bo». Wes improvement would 
enable n spectrometer +o still differentiate between valid pulse À and 


recoils without sacrificing all low energy valid pulses inthe process. It 





seems somewhat paradoxical that efforts should be made to extend the 
range to energies in the neighborhood of ł1.0 MeV from the upper cnco M 


? was aimed at extending the range og 


Maem thc original work of Batchelor 
his spectrometer beyond 1.0 MeV from the lower energies by eliminating 
the He recoils. 

A subsequent report {Persi will be written by one Health Physics 
Group of the Lawrence Radiation Laboratory, Berkeley concerning the 


actual functioning of the ?He neutron spectrometer, and departures (if 


any) from the predicted operational characteristics. 
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&. Previously constructed “Ee proport onal eot TE 


The design of the helium-5 proportional counter CUE a - 


a ker an extensive searcn of the technical literature. To aid Io Ae 


ጋ n 


ESMerciog or the counter tuoe parameters, 2 sunmecy of comperi ucc 


= 
۴۰۰0651651008 ر‎ filling pressures, and other pervaneny daca om aL “tle 
BEN r.10n&àl counters and lonization chambers consvrucucc Go Cave wes 
tabulated (Table Al). The Gata necessary to complete Table I is not 
ከ ቲቲ ር. Dil it 1S apparent that a variety of heavy monotonic gases aad 
Stabilizing gases nave been employed by experimenters. Counter dimensions 


and shapes have likewise been somewnat varied. 


Do lEvrence Radiation Laboratory usage requirements 
በስ - EM s کم‎ 

unc ue oroporcional counter will ultimavely be used for raga] 
Surveys by the Healtn Physics Devartment at the Lawrence Radiation Lavs 
EB The objective of this counter therefore, is vo be capable or 
E ine and accurately evaluating neutron spectra of energies to aDOUC 
14.0 MeV regardless of the neutron source direction. It is likewise 
necessary that the apparatus be capable of being readily transported and 
placed in relatively inaccessible locations. The usage requirements thus 
consist essentially of: (1) 4x resolution, (2) relatively small size, (3) 
the high voltage requirements remain within the capacities of existing 
Laboratory equipment, i.e. less than 6000 volts (an additional considera- 


tion that was imposed). 
gs Design variables 

Ia the design of a proportional counter, the designer has essentially 
five variables that are available to him, viz: shape (geometry), dimensions 
of the container (cathode), dimensions of the anode, gas composition and 
pressure, and applied voltage. It is to be noted tnat once a design ras 
5 è کے ہے‎ am AA * 7 2 157 oe ape 2 Ts E m . እ 1 ገ . ጫ 
been fixed, the only variable available to une experimenver 1S the applied 
voltage. Unfortunately, the five variables each effect tne tinal prop- 


erties of the counter and are essentially inter-related. 


e. Design parameters 
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TABLE A-1 
Physical Characteristics of 3Hc Proportional Counters and Ionization Chambers Reported in tho Technical Literature 
Type Counter Cathode Anodo Efícc- Opera- Gas Energy Fillins Gaa (atmospheres) at STP Experimenter, 
counter shape radiua radius tive ting multi- range He Stopping gas Stabilizing gae Date, 
(in. ) (in.) length voltage plica- (MeV) A: Total and 
in.) (volts) tion Qe) ል Kr, Xe. CO; Cli No  Preeeure Reference 
Prop, Cyl. 1,5 11.0 0-1.0 .071 1.45 1, 52 R. Batchelor 
(.0071) (1952) Ref. 1 
Prop. Cyl. 1,0625 , 0052 4.79 2600 ዴ10ዕ 0-1.2 0, 355 2.16 0237 2, 54 Batchelor, Aves, and 
(0.355) Skyrme(1955) Ref, 3 
Prop, Cyl. 1.965 . 0025 5,51 6000- 0-5.0 (1.0cc) 6-10 6-10 Bloom, Reilly, and 
8000 @STP Topple (1955) Ref, 2 
Prop. Cyl. 1,0625 ,005 4.79 4400 ~10 0255015355 4,46 ,0066 5, 84 Batchelor and 
(1.355) Morrison(1960) Ref, 4 
Prop. Cyl. 0.9375 ,002 9.00 0-17, 5 (0,976) 4.0 . 0418 5.02 A, R, Sayres 
(0,01335) 4.25 . 0445 4.32 (1960) Ref. 5 
(0, 992) 4,0 . 0418 5.03 
(0.0501) 4.25 . 0445 4.34 
Prop. Cyl. 0. 980 . 001 6,00 2450 4-5 7.0 (1.09) 4.50 . 0106 5.60 W.K. Brown 
(1962) Ref. 8, 9 
Prop, Cyl. 0, 480 .001 4-12 41200 (1-10) small 1-104 Mills, Caldwell, and 
a amount Morgan(1962) Ref, 25, 
Wang(1963) Ref. 16 
Prop. Cyl. 0.480 . 001 6.0 ~1900 (2) 4.0 small 6+ Texas Nuclear Corp. 
amount (1964) Ref. 26 
Prop. Cyl. 0.480 . 003 2100 (10) ? 10+ J. L. Friedes and R,E. 
0.480 . 003 920 (2) 7 2+ Chrien (1964) Ref, 11 
0.480 290) 1050 (6) ? 6+ 
0.3125 2006 5000 (9. 7) 93 10 
0.3125 ,005 5000 (9.7) 0.3 10 
Prop. Cyl 2.00 . 004 15.0 8,1 (1.0) 4.0 .01 5,01 Sayrcs and Coppola 


(1964) Ref. 10 
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TABLE A-1 (cont, ) 























Type Counter Cathode Anode  Eífec- Opera- Gas Energy Filhin; Gas (atmospneres) at STP Experimenter, 
counter shape radius radius tive ting multi- range He Stopping gas Stabilizing gas Date, 
(in. ) lin.) length voltago plicae (MeV) ————————— — —  TTsetal 2nd 
ün.) (volta) tion Que) A Kr Xe CO», CH, N2 Preasure Reference 
Prop. Cyl 0.4375 .0025 10.0 4000 450 $ 0.61 6.11 O. K. Harling 
(1955) Ref, 19 
0,4375 .0025 10.0 2500 11 ay 0,611 6.11 
d 0.4375 .0025 10.0 5000 93 D 1,00 10.02 
0.4375 .0025 10,0 4000 42 3 a 0.61 12.25 
0.4375 .0025 10.0 2400 32 Ü ቁ 0, 305 6,11 
5.93 DF) 
0.4375 .0025 10.0 750 1,0 Ü ir Pe 6,59 
nn 5.93 (5ኛን) ' 
0.4375 .0025 10.0. 10 12 Ü eds 6, 59 
5.93 (BF 4) 6.5 
0.4375 .0025 10.0 2000 2.4 ነ) تو‎ 9 
5.93 (BF 3) 6, 
0.4375 ,0025 10,0 3000 8 Gi a 59 
5.93 (5ኛ ን) 
5 : 4 6. 
0.4375 .0025 10,0 200 15 Gi nie 59 
0.4375 .0025 10.0 4500 390 ار‎ 0, 97 9,68 
0.4375 .001 6.0 1000 100 Ü aa ? ? ? 5.92 
jon EU 1, 375 «0. 14 9.25 6000 t. ” 3.8 0.07 ን ን Freeman and Weat 
: ' (1962) Ref. 7 
0, 46 : 
1on Sphero 3.94 0,295 N.A. 0.1-1.0 (0.46) 2.5 0.118 3.08 N. P, Claxkov 
s (1957) Rev. 13 
lon Sphere 0.511 .0197 N.A. 1600 0.2-0.8 2,.5-4 2,5-4 0,1-0,16 5, 1-8+ A. Jl. Abramov 


{1958) Rev, 12 





4 
m ma 


-90- 


I G እከ 11 v0 L2e Loree usage requirenmeavs, there e 
mechanical and electrical requirements that the counter must fulfi 


- 


Ecc O Cu Orens arei 


ደኀ ግ - Bd መወ , 7 ^ sa ce ር ውር ግ => re ጻጫሓዶ9ዊ ጫ ሥዛ ሥሙ 7 = 7 
۷۰۰۰٠٥۷۶2070 the iz requirement, the filling cas must ve capable of 


completely Scopping maximum energy product taxrticles within the Scc 


volume of the conteiner. 


(2) The well thickness of tne container must oe sufficiont to conten s 


Did ne weeSeS Under pressusc. 


(3) The gas fillings must not provide competing reactions that can be 


misteken for bona-fide events. Recombination snould be minimal. 





ሳጓ Gas multiplication to some derree is necessary to overcome electrical 
ic 


noise, and to obtain independence from the reaction site within counter 


(5) “Jitter-time" (the time it takes for an electron at a point Gaus 
cathode to migrate to the anode) must be reletively small, and within th 


resolving capabilities of the associated electronics. 


(6) Electric field distortion must be minimized, anà field grac ٦ 


equal distances from the anode should be ecual. 
(7) Electrical breakdown and "Corona discharge" must be eliminated. 
e. Effect of design variables on desig parameters 


The changing of any design variable has a profound effect 0 
design parameters. The interdependence of the parameters to several 
variables makes any proportional counter a compromise. The changing of 
ene wvyeariabues nas an eifect upon the paramevers as 7 


(1) Shape (geometry)--effects: lx response; field gradient distribution 


tt 


(thus "jitter time"); space charge distribution. 


(2.) Dimensions of the cathode--effects: amount, pressure, and composition 


of the filling gas; field gradients (thus "jitter time"). 


(5) Dimension of the anode--effects: field gradient ("jitter time” and 
gas amplification); corona discharge. 


(4) Gas composition and pressure--effects: stopping power; jitter time; 
gas amplification; competing nuclear events; relative cross ۃ٦‎ 


the composite gas. 





8 Selection of the appropriate desien variables 


25 The geometry selected was a cylindrical proportional coun Co 
cylinder waS selected instead of a sphere, due to she adverse effect of 

a concentration of positive space charces near the anode after the oceurr- 
E cof cas multiplication. While this condition lC Se ex M 
cylindrical proportional counter, the space charge effects are distributed 
over the length of the anoce wire and are not concentrated at a Sie Be 
point as in a spherical configuration. It was also considered thet the 
resolving time leg would be less (thus the counting rate higher) and 
multiplication more uniform in a cylindrical chamber. This selection, 
however, necessitated some sacrifice of the kx response of the proportional 
counter--a deficiency that was partially overcome by increasing the gas 
Eure. ihe field gradient distribution is radial along a central axis 
in @ cylindrical tube, except near the ends, instead of a direct radial 
GQistridution 2s in a spherical chamber. The end distortion of the field 
gradient of the cylinder can be eliminated, however, through the use of 


2 


appropriately designed field tubes. 


(2) The dimensions of the cathode were arbitrarily selected to be'a maxi- 
E Ge radius of 2.0 in. and an active length of 15.0 in. "These 
Gimensions were considered to be about the maximum that would meet both 

the portability requirement and still maintain a reasonable field gradient 
(thus jitter time) near the outer portions of the tube's radius. It is to 
be noted that the tube size, which has been dictated oy practical consicera- 
Prooumiikewise limits selections of the amounts, pressures, and compoSivion 


of the filling gases. The thickness of the cathode walls is determined 


Enc filling pressures. 
(3) Gas composition and filling pressures 
i.a.) The gas composition must be such thas: 


(5.8.1) The reaction particles from a 14.0 MeV incident neutron will be 
ScoOoped in a diameter of tüe cylinder. This will to a cervain measure 


کیٹ 


restore a portion of the hx geometry resolution. 


(3.6.2) The electron migration through the gas mixture should be suiiicienvly 





Shane 


Suificient to provice a pulse that is both proportional to the ence m 


ረ 


rapid that extremely long jitter-wimes are not encouncerec. 
(3.8.4) Competing nuclear events are minimized. 
i ኣ 


(2.5.5) A stabilizing gas be employed to absorb vhoton producecé from 


bremsstre inlunge during ges multiplication end to prevent spurious jee 


KJ 


(5.6.5) The sas multiplication (which is a Zunction of ihe press EDD 


ehe 1neident neutron and suiiicient to mask electrical noise. M 
تہ‎ ያ ። | un کی یر‎ 28,2 
MulGiplicavion Ot at least 10 is considered desirable, 7 ~ 2 
(3.0.) The selection of the gas composition was dictated by several 
TSc sO eS: 

(5.5.1) The cross section of the Ee (n,p)T reaction decreases by a factor | 


? 18 
or LO for neutrons at 0.1 MeV to LL MeV. It was decided tneat some 






compensation for this decrease in cross section snould be made if 17 
O1 higher energies are to be cetected in meaningiul numbers 6566 ከጠ 

reduction in efficiency at higher energies is likewise compounded M 
Lact that since the higner energy product particies will travel ss 


n the gas medium, their chances of intersecting a wall are in G 1 h 


W 


| 


and thus the probability of the production of a full energy count is 
4D ታ:; 


# 


Gecroascas was considered trat two 84 cmospheres of “He would compensate 


somewnat for tnis loss of efficiency at nigher energies. 


(3.0.2) Helium-3 is a gas with relatively low stopping power for charged 
particles. A heavy monotonic gas must be added to insure that the charced 
particles stop within the diameter of the cylincer. revious experinenters® 
have used argon, krypton, and xenon for this purpose. Argon ٦ | 
in this design for two reasons: First, the pressures oi argon Bee 1 
to stop a charged particle are greater than tnose of either Kr or Xe for 
ር ያ ዯ.ፒ.ር፣፣ ት Secondly, che Anya) reaction becomes Sigil .፡፡ዌግ 
above about 2.2 Mey. ^ ^79 Pulses due to this reaction can be removed by 
rise time discrimination, but the reaction data is incompleve Snom 
higher energies en argon (n,p) reaction occurs. ‘The disintegration cross ۱ 
section at a single energy of 14.0 MeV has been measured for argon, .n 
krypton, and xenon and found to greater for argon than either Kr and Xe. 

( The disintegration particles produced are not 700116767“ identific Sik 
reference 3L.) 


Xenon was likewise rejected due to the uncertainty of the type of 





Gisintegration particles producec, and abcenec of Catu conce n ደክ DE 


Qisinitogration cross section at intermediate energiocs. The rejocction 


f 


on 


xenon was Cciftficult since the heavier of the noble sases would ce ve 


بے 


meduined smaller filling pressures tran citner ar oTo To 0 


Krypton was tkus sclected to be tho "stopping gas". Unfortunately, 
ከ On likewise kas a D 8265576667 ناو 20030327 کے ہی‎ 
as an (n,o) reacsion by Shem 2 which produces a scattering of disintegra- 


۲۰۳۰۰۰۹76170760 energies over the entire energy distrioucion, Che e 
ZA ve 


E یت‎ : a 
lying in the 4-8 MeV region.” e S acvtempt will oe made to remove the 


ጫ 


ES INS produced fron krypton disintegration from the resultant spceutrus 


Memeo tne uncertainty of thcir energy dcisvriouvion. 


ENE Selected a heliun-j, -krypton mixture lt wes 96 60 . 

An maine the appropriate filling pressures. This wes accomplished py a 
۴۰۰۰۰۱۹ 6۶70< nethod of computation based upon detvcimining uno lengen 
EE ck of a proton navins the energy of the average cosine of Scavuverins 
angle (Lab System)from a 14.0 MeV incident neutron. When the length of 
pos as equal to the radius of the cylinder, it was considered thay une 
pressure requirement had been met, and that a configuration trat coula 
approach the desired 42 resolution had.been accomplisned. It was found 


E 


that 10 atmospneres of krypton plus 2 atmospheres oz "Le meu this 


"ን ን 10:0. I 
(3.0.3) A stabilizing gas was considered essential and 0.5% CO, was 
K fOr this purpose. Carbon dioxide was selected in 60 ٣ 
v 
E : - 


which has been successfully utilized oy Batchelor, Freeman and West, 
and a in their most recent counvers. Uno reason to kihisi To nS 
is that our counter is designed to evaluate neutron spectra to 14.0 MeV 
(considerably above the energies evaluated by previous experimenters), 
ano it was considered that the introduction of a source of additional 
recoil protons (from the methane) into the tube would provide unnecessary 
complications--especially in view of the additional proton background due 
to the krypton disintegrations (which wes not a problem at lower neutron 
energies). 

(3.b.4) For comparison purposes, the cross-sections of tne various gas 
reactions are plotted in Fig. Al as a function of neutron energy. ‘ine 


ESS section for the CO, was Ootained by adding the energy dependent 
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(barns) 









c > | 
ር መ ር He (n,p) T j 
O I^ (6) ?He non —elastic 
© 
ርን 
uo 
ሀን 
O 
O O| 
ae 3 
(7) 9 (n. d) He(n,d) D / 
بت‎ a) Kr (n,p) Br 
(oS 
کت‎ 
= | 5 He(n,p,2n)!H 
| | 
| 
0.001 ' 
Q.1 0 | O | OO 


Neutron energy (MeV) 


ለላህ. 37191 


Tic TAI. Neutron e cert osa ccm 


Source OF Data 
(Ime @aleulated from 6 2) Garo TOL OS vac (Rei) 
TOE ወ 4 E 
(3-7) Ref. 18 
(8) Ref. 32 
(9) Ref. 18 








100 ; | 


| | | 


| | | 

| 

' 
۱ | 

(ጋ | 2 
| | (5) c Hex2 
| = 
| (4) v elastic 9 He (n, n)? He xe 
3 : 

1 É Ge a p) PIT x2 
| l 
| ህክ - መ: 4 ። 
۱ ; i [| | (6) "He non—elastic x2 
| | 88 IL 






















== F 
C 


(barns) 





| ne, CO, x0.05x12 V Y 

= ` م‎ vY h ^ سس‎ 

6 | | "eA i 

= ae 

A lovee Ate "m 

ርን ዕ | (n.p) r n.p rx a 

ሃን "m | 

o 

o | 

- f 

© | ; | | 

(ሀገሮ Š He (n,d) 0x2 

> (n,d) | KOC on) 

=  SHe(n,p,2n) H x2 

| نہ 

መ | 
| C 0.01 | 


| 
0.001 ot | x | 


on os 8 10.0 100 


Neutron energy (MeV) 


ለላህ - 37292 


Hig. AZ. Relative Neutron Reaction Cross Sections in a 
Gas Mixture of 10 Atmos Kr. 2 Atmos te, 
and 0. 5% CO, 


> s/h 


. 
-—— — 


- 








toča. Cross section of carbon to twice the total cross section of oxygen. 


in addition, a relative plot of the cross-sections is presented to o G BE 


tae collision provabilities "seen" 


SC. - ሠ= መጣ ہہ سر ہے ہے ہہ ہل کے‎ ean > 
py S 100-0 00001 AS lC CGO ce Cage 


(Fig. A2). It should be noted that the disintegration cross section or 


f) 


aoe UE ግ -- : 5 - ሥጓ =.ጫ = ሥማ ~ ግ e `... f^ መ= yA S ወሙ d al‏ سے መ‏ ^ سے 
mW CON 2c LO aumospheres prossure boconos & factor of a GO‏ 1 
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5yüeroclod qu nicnor energies, anc. probably will aevcrine vane ک‎ 


prr m Ebo ያሬ መመላ ገር በር On. Us S typo coungen, 


(4.0.5) Tne dimensions of the anode were now considered and several appli 
] 


| 


voltaze and Wire size combinations vere attempted. Jitter time was 


LH 


calculated by use of a "weighted" electron mobility factor (preseateG on 

page 19 ) and found to vary from 10-22 useconds for applicd voltages of 

6000-5000 volts. The gas multiplication was calculated from extra B 
0 


LOU 
of Brown's Cata™ and found to be approximately 10-15 for tle ao ا‎ 


ጣጫ رب‎ 


of .voLtases, using a j mil diameter wire. The jitter time is con ns D E 


[ 
longer than desired, but increasing the wire size (for fixed cathode 


not appreciably decrease the jitter time in ihe voltage ranges Com as 
The jitter time is relatively insensitive to changes in 26وہ‎ 1 
therefore if we seek to reduce the jitver time, che only Leas BL a 
“4S to reduce the diameter of tne counter and increase the gas ٦ 
pressure. This alternative is undesirable since an increase i eee 
pressure would increase the relative disintegration cross section om M 


krypton--a problem that is already larger than negligible 0و‎ 0 ٣ 


(2.5.6) We must therefore accept the fact that the counter will be i 


sonceresrvely slow, and design uoce cleccvoutes 2ccer 1 ae 
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S ^e, m E مل‎ `Y 2 ብልጡ - -“ መኣ بس ا‎ 
fn Of the Associated x. ecuroni G Go ህክ 





s e ren 
a. Design Parameters 
ena s. = ከፕ መ 


Tne electronic components of the “He proporciona ro ٦+ EM 


s I‏ سا 
S‏ تد ተመን...‏ 1 


- 


5 ሥነ 


ነ የሮ! ዝቃ, meet several requirements. The majorivy o2 
Memeo are common to proportional counters in general, However. 6 - 
EE J] Unique problems associated with the broad range or 77٦ 
of the pulses (essentially 3-22 usec is our region of interest). The 


@esigem parameters include: ; 


(1) The pulse produced by an ionizing event must be linearly amplificd 


CS Capable of accurately portrayins the energy of the ionizing event: 
(2) The electronics must have a low simal to noise ratio. 


(5) The high voltage supply must be well stabilized, since variations 
BN De hign voltage will cause a direct variation in tne pulse height. 


^ . 


(4) The electronics must be capable of distinguishing between a desired 


event and an undesired event (recoil). 


(5) The electronics must be sensitive to the relatively lerge range of 
Mase cimes, and to accomodate pulses that vary in botn energy and rise 


time. 


6) Provision must be made to vezmit variation in: the high voltage 
ae ° 
EE the counter field tube settings, the rise time discriminacor 


ጫ 9 (0 "5 


Es. and the time constants of the double Giilerentieting nevwork. 
EN "ecuvtronic design 


(1) A block diagram of the electronics arrangement constructed to meet 
the above parameters is indicated below. The basic arrangement of 
Pemponents is similar to thav described by Sayres and ےت‎ with tae 
exception of a new Peak Sampler Unit and a Time to Height Converver Unit 
which were designed and constructed by M>. R. X. Brown of tne LAL Physics 
Tecrunical Support Group. The purpose of the Peak Sampler Unit 1s to 
provide a gating to the Multichannel Analyzer at a time ween ihe puise 
has reached its maximum value. 


ጣሣ ° 


(2) To explain the operation of the electronics, consider Loe forme clon. Cam 


pulse in.. the proportional counter as gne result oil a coll L RK m. 


ESS Th cone ol tie contained sas molecules arca سا‎ co, )- Tae 
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Fig. A3. Block Diagram of Electronic Components for Rim a 
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ጫ ." x e 7 ሓ > کت‎ e es = * መዓ m - 2 2 ግ : M و‎ I. ጫ S بے‎ ፍጫ Jg መፍ a E 2 UI 
collision may resulv in either a Cesired event, 1. e. G GS c 0 


መ: ሠና“ E >¬ - e 2) ام تے ا مت‎ "T Res کا ند‎ ^ e E 
boe “He nucleus, or an uncesirced evonií--b5e ጋመ 1፡7 R ጸ 


ወ 


r = > A جم مہ سا‎ ü -< ` “ታነ ےم‎ -} ye ve Vs x em-e ሥነ کٹ‎ - "ግታ ግ \ 
19:2.) frhe electrons that wore libarated as tic ona doc USER 
+ ہہ ہہ‎ c: Š J- mm cas arce el e AR ٦ تہ‎ AA ہے ہل‎ መጋር. 21 ہے‎ ^-^ ER m am da °" m m 
DONS S UC Sas i eC alca DOLI 8 28 ٤٤ مہات سی‎ 1۷٦ ana a 


ہا کے 


@herse is ceposived on the anode. .The size or the chr ciis denendene 
nen che energy or the particles that are stocned Wltnia thew 8 . 
s S amount of ges multivlication in the area surrounding peewee. 
(For ges multiplication above about l0, the exergy spectrum is undistorved 
Ene position of the ionizing event within tne ube"? providing, ٦ 
EX On OL positive ions in tne sneata surrounding wae ancc] r R a 


been achieved.) The electrical pulse bevins to accumulate when the first 


© 


۳۰٠۰٠٠٠ reaches tne anode, and reaches 3 maximum Ween the OS 1-۳٦ 


۲٣۰۰۰۰٠۰ ۰77۷65 at the center wire. Depending upon one re O even. 


۳۰٠۰٠۰٠٦۰. energy of the incoming neutron, tne radial distance 170m) 
EN SC resuiving vraex orienu&vions Witch respecs vo whe anoce, 7 
Charge build-up period may vary froma fraction of a micro second to as 
ከ222 wseconds for a given voltage gradient Gistributvion Within che 


ጫ 


" ከ... Tre actions of the following component units begin 06 


09 


5 በሮክ ሥ2ሮ begins to appear on the anoce, anû continue Until 60 


Ec hes oeen collecte 


(2.b.) The incoming pulse is amplified through a charge sensitive 


Preamplifier and an Amplifier. 


(2.c.) Tne pulse is sent into two parallel channels after leaving the 
ee. One portion of the still rising pulse ሠ. ድ ር ይክ L 
the Peak Sampler Unit where it accumulates to @ maximum--at which cine 


it is either "gated" into the Multichannel Analyzer, or is discarded. 


feo.) The other portion of the pulse is sent to a Double Ditterenvieting 
Circuit (figure A-4.)- where an output pulse is obtained whose duration is 
dependent upon the rise time or the input pulse. (As ean approximation, 


the output pulse from this unit represent tre “rate oi cime” since 


Semele RC circuits are not true cifferentiavors. Tne Gime OCCT 
line crossings of the doubly differentiated pulse is thus relaved to ‘tne 
rise E The exact relationship between the Gifierentiating time 


- 


Gonstents end the rise time is not critical since we are only interes vec 


Meesctcing a discriminetor value waich eliminates -eall oulses having ea 
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መ cee ግ ہے‎ a 2ب یں‎ NC E : ዘ ወ MEE ed 
ናፕ ከ ቲ ክር... ..1...1 ርፌ ር. (ርመ uC ii Nec Vous Cold Cs 257465. POSS WES EC 
سط م‎ er EU dae 5 r. a s MOT 

ter than tnis value. ae SILI OOM ope e S208 Jg‏ 7ج ت5 ۱۷۳065 تح ا2ت 


by eee. 2 ٦ - ን ሚፍ E መመ ታን 5 ር A ዋሃ” ገ — ہر‎ 
c 0601760607 06ک‎ ubtL.l2.36 a sical cCenereavos vita See 


known wise times. The shape of ihe input anā output pulses are shoma on 


وس 





0 
gig. A-4. Tne differentiating time constants arc cqual to eacn ou R 
each of the three settings (1, 2, and 5 usec), and can be varied depending 
upon the desired output. The counter normally operates with tne cii EEE 
tiating time constants sev at 2.0 usec. For an analysis or 002200 የጸደይ 
spectra, the differentiating time constants are 5601 aco 1... ዚኃፎይር፡::2.. 8 
AP d SS P t > 2 
che higner 5.0 usec time constant when we are primarily 3126 2276606 ٣ 


ኣል] 


ne plebex end oi whe رت‎ 6 77 50606073۰ 


cr 


(2.e.) The output of the Double Differentiating Circuit is 2۶ل‎ 
unit called a "Pulse Height Discriminator". Tne purvose of tnish r 
to convert the input pulse to a square wave, and from this وہ‎ MZ 
generate a "start pulse" anda ። pulse" as the input double differen- 
tiated pulse crosses tne zero potential axis. The electronic arcc 
ror this component is snown in Fig. A-5 along witn the circuitry Doce 5 የየ 


COE DIOGO bnew VOL Dies source: 


(2.£.) Tae "start pulse" and the "stop. pulse" are then directed to & 
ime-to-Heignt-Converter unit. A block discription and a 1۸8 99 
diagram is displayed in Fig. A-6. The "start pulse" sets the Bistable 
Multivibrator thus initiating the ramp generator (consisting of a constang 
‘current source charging a capacitor). The charging process continues | 


until it is stopped by the action of the "stop pulse”. The "stop pulse" 


is red into a Monostaole Multividrator which produces 4 2.5 1 
oulse. "This pulse turns off the constant current source, ይስዕር 


vermminates tne charging process. The heignt or the raro at UvDeN NEM 

waere it is terminated determines the height of the 2.5 usec mi dci TEE 
pulse. due ሽሽ... is rapidly discharged, and the Start Muitivibrator 
1S resev ror the next pulse. The output of the Tine-to-Heignt Cony. ማየ 
ts 5005 8a pulse of تک‎ such tat it is directly related to i sS 


it output--thus directis 


H 


Goss Cr the double Gir crentiacin= ec: 


c 


related to the pulse rise time. The pulse width of 2.5 usec So 7 


Significance other than a convenient width for the ۶0110112-8 0 


4 


A circuit diagram of the Time-to-Heignt Converter is snown in Fig. A-~-f. 
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Fig. A6. Block Diagram of Time to Height Converter 
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FIG A? ELECTRONIC CIRCUITRY, TIME TO HEIGHT CONVERTER 
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(2.g.) The output of the Time-to-Heignt Converver is then sent to a 


- 


Single Channel Analyzer (Fig. A-&). The Sinele Chennel Analyzer can se 
set to accept all pulses above 2 certain height (cine) or to aceon: 


— 


puises on a winGow basis. The analyzer settings are adjustaole. xunbers 


Y 
+` 


s i. A-8 refer to standard IRL årawings . Tmerr o nm. 


EE ህጋር፦ produces an output pulse only if the input emolivude 25522. 


Mer winaow Setting. 


(2.h.) The output pulses of the Single Channel Analyzer (rise times 
Within desired ranges) &re then passed to tke trigger input of the Peak 
Sampler Unit (Fig. A-9). The original pulse has been building up to a 
maximum and when the "trigger simal" arrives from the Single Channel 
Pee yzer, the pulse is gated into the Multi-Crennel Analyzer. If vne 
trigger signal" does not enter the Peak Sampler Unit, the pulse is iod 


poscround. 


* Data on the units of Fig. A-8 can be found in neference 25. 
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Fig. A8. Block Diagram, Single Channel Analyzer Unit 
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AN = ^T wm - ne ቁ کہ م‎ 
Avpencix B. h CET ECO melis 








Consider voe General CE Ones Porvicle of macs mane Vy +٦ 
whe 
Wiun & particle of mass 2... assumed Xo Dc et r6sv. ጋ ከ100. m 6 5 
6 . . ۸۲0667۰-7167 بک‎ ۹۶۶5 ene m, with veloc a 


and. ME leave the reaction site av angles 5 end “respec aa 


+ 4 
- سل 


L5 


s ግ 


Utilizing non-reletivistic dynemics, tne laboratory system or COC me ስከ am 
and mass numbers of masses, We can write tne tnree general equation sS 
tae conservation or Monencum and euere,. 


ጫ 


ር ን 118.7... ገገ 12. 


I. V. 2 nm. V. COS -m v, cos V 
mE S ly ስ. በክ ( 


.. C pare ኒ' ር 
0 = m, V, Sin ò- mV, sin Y, (B 2) 


S 


b. For wae Conservation of 7 


(5-2) 


where Q is the energy equivalent corresponding to the mass difference (in 
non-relativistic particles) between the initial masses (before collision) 


and the resulted masses (after collision). i.e. 
Q = (m. + m) - n + m, ) (5-5) 


Note that Q may be positive or negative. Q values are tadulaved 


TA 


ማ ው . A / 
5701676 readily መ 1 1220.ሀ/:7ን፡.ፕፒ ያለቃ ۹۹پ‎ ٦ . <B-4) 


Equations (B-1), (B-2), and (3-3) may also be written in terms of 


rey by making ہہ رتو وک‎ aN ee 
energy by making use of the relationship BD = 5 mv" of mv = 1 
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ሚመ uU E eM . کے‎ መመ 
ersero eo ace sae 


constant {2 factor in Eos. (B-1) a: 








j E ] 1 
Y^ TS - ^4 ማሜ c ~a Tn = L 3 \ 
ደር ከ siu. 1. COS yy (2051 52 
EE c ncs n E لس‎ e 
O = | yn T = cm E [m - v 
— ءل وی‎ Sus V, s=) 4. ! ے١‎ DLN n m (5-2 ፡፡ 


In the following analysis, we small work with energies rather than 
psereles. Noto tuat by the sixultancous Solutions or vto QUY s 


equations, we may eliminate any two raremeters. 


Cace D. oin gue Doubs NE 


ll 
ri 


Meo us irst consider the case of Simple scavtering 7٦7 


5 
and m, = m,- It is readily seen from Eq. (B-L) that in this case sqr see 
Tv 1S desirable to derive an expression for the energy imparved to the 

Stavlonary target na 
the above mass substitutions into Eqs...(3-l.a.) and (B-2.a.) transposing 


፥ <= DANI 2. 3 a S “ሠ. کے ہہ‎ 
=n terms of the incidenti particle cnc 099 ame 


‘the last term of Eas. (B-l.a.) and cither term in (B-2.a.), squaring 


Iibeva equacions and adding we can eliminate the angle 5, to obtain 





j T ኒታ መመን = JS 
በ - Tn IO Ped COS m 4 = M, š 
2 T ze ፡ 5 ጋ. ተጨ 


Solving Eq. (B-3.a.) for ጅጋ and subsclouviHs cals + G n Q ንዓ 


above relationship, one obtains: 


923 





> fsd 
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۰3۰ n Sides anc solving for =), Vue vase 
በ 
Sant 4 ጊ 2 m a CJ ~ — 
T - < بے‎ 3 ያ COS ٠ (5-5) 
ہے‎ 77 + “L = RE سرچ‎ 
ርር ار‎ 


Therefore, the energy of the recoil particle is dependent upon (a) 


مر 


the masses the incident particle and target particles, (5) the energy oz 


ርን 
ct 
D 
(D 
cl 
0 
ትን 
C) 
(3 
ci 
Ij 
m) 
! 
ረ፤ 
b 
Q 
፦“ 
(D 


A M QU CT REM S IUE‏ ]1 اگ دے D‏ 203 .` دی 
paruleleswasd c) P OONSCRUS ODE‏ 5 


= 


rakes to cke direcion oi Cre Tinei den. o Lele: OTE ximum 2 2 ٦ 


۰ 


"s M 
Ea er 


(n 


OOVLOUS 21፡2 ሮጋ፣ ጋዳ ከ ን ር 


ኑን 


en = 0] 853.1. impart the 27 
ELECT TO COC CCI GICLEE 

Our counter is designed to analyze xeutvrons of energies aos 
14.0 MeV. ? is of interest to devermine one maximum energy TA ٣ 
imparted to a potential recoil particle. in tais case, ler Oco ú 
ipn E = O. Using the atomic masses of tne neuvron and the Led ስርርር 


nucleus, and tne chemical avomic welgnts for kryoton, caroon, a as 


respectively, we can substitute these values into Eq. (3-5) to obtain: 


Maximum Recoil Eaergy 


General for 14.0 MeV neutron 
dT 8 ue recoil cC hors l0.5 MeV 
max 2 
Bo UE carbon recoil 0.284 E. 3.95 MeV 


ë J i É : M 
s Bur oxygen recoil 0.222 By 23.31٠٥ 


7 ° I TA ے2‎ - 
L ZNN CT TUCO eco) 1 0.0400 سک‎ 0.05 Mew 


Case II - "Ee (n,p)T reaction 





2 


In the case of the “Ee (n,p)T reaction we can write the equations 


for conservation of momentum and energy as follows: Let us idenvily m, 


aod E, as a neutron having mass of unity and energy M mass m, 12:008 | 


ud 
po tne “he molecule having s mess ot 5 


" 


Snall 1006 m es uhe 


E 


triton having a mass of 5 and energy Dns Inc ር | وت‎ à mess of 
Meloy ክን ነኸ ፐን. my and | Mave "an eaccey 4 Makirsg these sSuUuDE M 


wLons into mds. (B-l.a.), (S al) and کہ‎ yields: 


. 


TEE — ara tagen tit کے و‎ 
ee = vi I. COS > ጭ /ጄ. COS Y E (3-1 8 ° ) 
N ፌዛ ፍ i L V L (oe 
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da E K ጣ “ ኻዬ مب‎ Ld . f - 
7 = መጨ - - ውና ቁ Ai. 7 کے ہہ‎ " nee -- E E i" ` - . 
ES DOO. in terest.an the ہ٥ بت‎ 62720302802001 ox rico Sie sS O M 


Gave an expression for the angle C, between vae proton vector E. anà the 
A 

- 7 ሬ - መመ ! jos 5 ግ መጻ “መመ = سم‎ p “ል کے ہج‎ ጫነ - = 
Seton vector Mme Thais expression cen be ostained by squering Ecs. (3-i.b.) 
and (B-2.b.) and adding the two equations. This yields: 

; pL ንን o i sŠ = sinV 

7 2 V9 ES Es (cos V cost inV siné,) > 

Recomizinz thas the bracketed tarm ና : B ኣ/ < 
RECOZNIZINg wnat tae bracketed term is equal to: cos ( ረ o, ) = 


cos Q,  , we can solve for cos OL, 





ማሣ 


The solution of equation (B-3.b.) for Zand the substitution into 
i | 


۳۶۰۰-۶6086510 Lor cos œ, yields: 
+ 25 
WO T 8 
COS Q, = ( — ) (5-5) 
” ራታ ኤቢ Ss 
{9 By Bp 


Note that œ is always obtuse since Q. D is- -a positive o niin 
2 


5 . Sh M 
a | Q TO MeV 


- 


mile Sin C will also oe of inverest in Our CORDLESS 6 ብ - 


easily derived using the Pythagorean وت‎ to produce 


(ia (2 ay -F 


SNC —— MÀ (3-7) 


Ext E 





بط جب 


5 
Case III ^He(n,8)D Reectio: 


cdueteron producing reaction, lev us 2Gentaly Gae vcore ን ከ ከጌ. ጌጉ ሱሱ‏ ے۶۰۰ 


"> ክርን 1125: Ey aad unity are egein subs civuvec 


e ۰۱٠١۱۷ Of mass 2; P By, ê anc mass A SS c = 3 


። ያ/ 
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nuance are Uti lized for , and m,- Suossivution oz 
"T ^r 
norosos. [5-l.a.), f5e2.a), end (8-5.8.) yieids: 
ከ መ a : Eo a etu 
US ہم نا‎ cos à + 92 Boy COS Ye 
ያ 5 (UJ 
ህን እ irae رر ہت‎ ٣٦ Saas 
ኒ Dx. 0 D2 ፦ 
E ` = 5 4 E 
N "5,0 Dl D2 
We are arain intercsced in tae መ ።.ሮ c b° 2 i 
ገ ከ. OCO 2۰513 On pz ۱ت‎ 8 8 +6 CN 
E -2 5 A ብ. 
N DL D2 
cos 2; = — 
Lig E 
d Dl. “D2 


Solving (B-3.c.) for 5. 


+ EQ, + 5. 
(Qa EDN oi) 


ከ T 7 
P OE 


COS Q = 
ፌ 


ከው ታእ ውም ን ን ንያ. لا‎ p E O =F Ç which is negative 
7 D2 mod ( p ); 
m ሥ ° 1. mA +m m 
ጄ E © C. 256 acute; ana 1i E FO Q 
DENN E EE QE : 32 * 955 7 


- 5.27 MeV and is threshold &t h.56 MeV 


can be determined as in case II 





cae two deuveron 


ጨር SUDSTLUUTLON into tne value ٦ c. 


then COMES 
L 


፲, is obtuse.‏ ر۔ 
Qc ۶۲‏ 


“>= ) - i 
CLES Vol 


(B-L.c.) 


(3-2 P 


(3-5.c.)] 


HW 


= 


| 


(B-9 








y M ty‏ س سو 


እ ን eee 
4 ^ لو اد کے‎ È 


ማላ መመ መዛ TET دا‎ E — ہا ہے‎ መቹ ያ ቸው / ጫ ለ መው ا‎ r ባሟሙጭሎ - ል ብቅ 1 መጨመ ہے‎ ው / : - پ‎ 
VIDSINLDULLON OF ےراھمت‎ AND ANGULAR DEPZNDENCH TN Or; Z K ንኑ ህህ 


aay S کے‎ ደ 24. መተው m , Pod d an = ipm s - ے١7‎ = ff Ves uu گے حم‎ 
A. Reason for tne Variation of Cos © مم عے۔‎ TC Os b a D 


A ac A ሥሪ ሔ.ሴ 


ሬት 1 - m- LE ሓዛ m E. ہے ہے‎ e san wien oy ^ 
NSEC ME EOD LT os IO 


^ . - 


Zn spaerical coorciantes, 5. 5152052205 02. እረ جج عو‎ ን 





(C-1) 


li we wish ቲዕ cover every Dossin -C 67 7)۹ ۶٠٦ oe sphere; 
LN see that the differential erea is not composed of merely do ard 
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This report was prepared as an account of Government 
sponsored work. Neither the United States, nor the Com- 


mission, nor any person acting on behalf of the Commission: 


A. Makes any warranty or representation, expressed or 
implied, with respect to the accuracy, completeness, 
or usefulness of the information contained in this 
report, or that the use of any information, appa- 
ratus, method, or process disclosed in this report 
may not infringe privately owned rights; or 


B. Assumes any liabilities with respect to the use of, 
or for damages resulting from the use of any infor- 
mation, apparatus, method, or process disclosed in 
this report. 


As used in the above, "person acting on behalf of the 
Commission” includes any employee or contractor of the Com- 
mission, or employee of such contractor, to the extent that 
such employee or contractor of the Commission, or employee 
of such contractor prepares, disseminates, or provides access 
to, any information pursuant to his employment or contract 
with the Commission, or his employment with such contractor. 
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